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ON THE STELLAR POPULATIONS AND EVOLUTION OF STAR-FORMING GALAXIES AT 6.3 < z < 8.6 
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ABSTRACT 

Observations of very distant galaxies probe both the formation and evolution of galaxies, and also the nature 
of the sources responsible for reionizing the intergalactic medium (IGM). Here, we study the physical charac- 
teristics of galaxies at 6.3 < z < 8.6, selected from deep near-infrared imaging with the Wide Field Camera 3 
(WFC3) on board the Hubble Space Telescope. We investigate the rest-frame ultraviolet (UV) colors, stellar 
masses, ages, metallicities and dust extinction of this galaxy sample. Accounting for the photometric scatter 
using simulations, galaxies at z ~ 7 have bluer UV colors compared to typical local starburst galaxies at > 
4 a confidence. Although the blue colors of galaxies at these redshifts necessitate young ages (<100 Myr), 
low or zero dust attenuation, and low metallicities, these are explicable by normal (albeit unreddened) stellar 
populations, with no evidence for near-zero metallicities and/or top-heavy initial mass functions. Most of these 
galaxies are undetected in deep Spitzer IRAC imaging. However, the age of the Universe at these redshifts 
limits the amount of stellar mass in late-type populations, and the WFC3 photometry implies galaxy stellar 
masses ~ 10 8 - 10 9 M Q for Salpeter initial mass functions to a limiting magnitude of M1500 ~ —18. The masses 
of "characteristic" (L*) z > 7 galaxies are smaller than those of L* Lyman break galaxies (LBGs) at lower 
redshifts, and are comparable to less evolved galaxies selected on the basis of their Lyman a emission at 3 < z 
< 6, implying that the 6.3 < z < 8.6 galaxies are the progenitors of more evolved galaxies at lower redshifts. 
We estimate that Lya emission is able to contribute to the observed WFC3 colors of galaxies at these redshifts, 
with an estimated typical line flux of ps 10~ 18 erg s" 1 cm" 2 , roughly a factor of four below currently planned 
surveys. The integrated UV specific luminosity for the detected galaxies at z ~ 7 and z ~ 8 is within factors 
of a few of that required to reionize the IGM assuming low clumping factors, even with no correction for lu- 
minosity incompleteness. This implies that in order to reionize the Universe galaxies at these redshifts have a 
high (~ 50%) escape fraction of Lyman continuum photons, possibly substantiated by the very blue colors of 
this population. 

Subject headings: early universe — galaxies: evolution — galaxies: formation — galaxies: high-redshift — 
intergalactic medium — ultraviolet: galaxies 



1. INTRODUCTION 

Galaxies evolve very strongly as one observes them during 
earlier cosmic epochs. At some point in the past we should 
begin to witness the periods during which galaxies are very 
young, likely having formed no more than a few generations 
of stars. Strong evolution is observed in the integrated galaxy 
star-formation rate (SFR) density, which increases by roughly 
an order of magnitude over the redshift range z ~ to 1 , and 
reaches a peak at some point between 1 .5 < z < 3 (e.g., Lilly 
et al. 1996; Steidel et al. 1999; Dickinson et al. 2003; Hop- 
kins 2004; Reddy & Steidel 2009). There is a marked de- 
cline in the ultraviolet (UV) luminosity density from z ~ 3 
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to z ~ 6, implying a lower SFR density when the Universe 
was < 10% its current age (e.g., Sawicki & Thompson 2006; 
Bouwens et al. 2007; Reddy et al. 2008; Reddy & Steidel 
2009; Bouwens et al. 2009a). 

Star-forming galaxies are readily identified at z > 2 as the 
redshifted Lyman-break (and the increasing opacity in the 
Lya forest) moves through the optical passbands for 2 < z < 6 
(e.g., Steidel & Hamilton 1993). The broad-band, multi- 
wavelength photometry of typical ("L*") Lyman break galax- 
ies (LBGs) at z ~ 2 - 3 shows that they have rest-frame UV 
colors consistent with local starburst galaxies, whose light 
is dominated by young stellar populations (on the order of 
10 Myr to 1 Gyr) with modest amounts of dust attenuation 
(e.g., Meurer et al. 1999; Papovich et al. 2001; Shapley et al. 
2001, 2005; Erb et al. 2006b; Reddy & Steidel 2004; Reddy 
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et al. 2005, 2006, 2008; Reddy & Steidel 2009; Bouwens et al. 
2009a, and references therein). 

Star-forming galaxies have implied stellar masses of 
10 10 M at z ~ 2-6 (e.g., Sawicki & Yee 1998; Papovich et al. 
2001; Shapley et al. 2001, 2005; Yan et al. 2005, 2006; Erb 
et al. 2006a; Fontana et al. 2006; Reddy et al. 2006; Overzier 
et al. 2009), with some LBGs having stellar masses as high as 
- 1O U M (e.g., Shapley et al. 2005; Huang et al. 2007). At 
higher redshifts, 4 < z < 6, the LBG population shifts to bluer 
rest-frame UV colors (e.g., Papovich et al. 2004; Overzier 
et al. 2009; Bouwens et al. 2009a). While some LBGs at 
these higher redshifts have implied stellar masses as high as 
~ 10" Mq (e.g., Yan et al. 2005; Eyles et al. 2005; Yan et al. 
2006; Eyles et al. 2007; Stark et al. 2009), the stellar masses 
of L* galaxies decline with redshift to values of several times 
10 9 M Q at z - 6-7 (e.g., Stark et al. 2009; Gonzalez et al. 
2009; Labbe et al. 2009b). 

While our knowledge of the properties of galaxies at 2 < 
z < 6 has grown substantially, at z > 6.5 the combined Lya 
and Lyman continuum break moves beyond ~ 1 fim, into the 
near-infrared, where the terrestrial background limits astro- 
nomical surveys to only the brightest objects. While several 
recent ground-based surveys have searched for z > 7 can- 
didate galaxies using deep T-band imaging at 1.0 /.tin (e.g., 
Castellano et al. 2009; Ouchi et al. 2009), these have pro- 
duced few high redshift objects. Even the deepest NICMOS 
exposures are limited in sample size due to the relatively large 
point-spread-function (Bouwens et al. 2008; Gonzalez et al. 
2009). 

The recent installment of the Wide Field Camera 3 (WFC3) 
on the Hubble Space Telescope provides the ability to ob- 
tain deep images in broad-bands from 1-2 fim. Recent 
deep WFC3 imaging in the Hubble Ultra Deep Field (HUDF; 
Beckwith et al. 2006) allows for the identification and study 
of galaxies selected primarily as LBGs at z ~ 7-9. Sev- 
eral groups have capitalized on these new WFC3 images, and 
have studied the luminosity densities, average stellar masses, 
morphologies, and inferred ionizing properties of galaxies at 
6 < z < 10 (e.g., Bouwens et al. 2009c; Oesch et al. 2009; 
McLure et al. 2009; Bunker et al. 2009; Yan et al. 2009). 

Here, we use a new sample of galaxies at 6.3 < z < 8.6 
selected from an updated reduction of the WFC3 imaging 
combined with existing imaging from 0.4-1 from the Ad- 
vanced Camera for Surveys (ACS) on board HST and deep 
Spitzer Space Telescope observations (3.6-4.5 pm). In §2 we 
describe the data and galaxy sample selection. In §3 we com- 
pare our sample to those obtained in previous studies of these 
WFC3 data. In §4 we investigate the UV colors of these ob- 
jects. In §5 we compare the full spectral energy distributions 
(SEDs) of these objects to stellar population models, putting 
robust constraints on their stellar masses. In §6 we discuss the 
effect of Lya emission on the galaxies' SEDs of our sample, 
and in §7 we discuss the implications of this sample on cos- 
mic reionization. In §8 we present our conclusions. We de- 
note photometric magnitudes measured in the ACS F435W, 
F606W, F775W, and F850LP filters as B435, V 606 , i 175 , Zssp, 
respectively. Similarly, we denote magnitudes measured in 
the WFC3 F105W, F125W, and F160W filters as Y W5 , /125, 
and //160, respectively. Where applicable, we assume a con- 
cordance cosmology with H D = 70 km s" 1 Mpc" 1 , O m = 0.3 
and fl\ = 0.7. All magnitudes are reported in the AB system, 
where m AB = 31.4-2.5^(^/1 nJy) (Oke & Gunn 1983). 
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FIG. 1. — Photometric redshift distribution of the high-redshift sample. 
All galaxies in this sample have S/N(J l2 s) > 3.5 and SfN(H 160 ) > 3.5 and 
Vi > 0.6, where "Pg is the photometric redshift probability distribution func- 
tion integrated from z = 6 to 11. The black, blue and red histograms show the 
distribution of best-fit photometric redshifts for all 3 1 objects, the 23 objects 
in the z ~ 7 sample and the 8 objects in the z. ~ 8 sample, respectively. The 
blue and red curves are the summed photometric redshift probability distri- 
bution function of all galaxies in the z ~ 7 and z ~ 8 samples, respectively. 

2. DATA AND SAMPLE SELECTION 
2.1. WFC3 Observations and Reduction 

The HUDF was observed by WFC3 from 2009 August 
26 to 2009 September 6 as part of program ID 11563 (PI: 
G. Illingworth) in three broad-band filters, F105W (16 orbits), 
F125W (16 orbits), and F160W (25 orbits), with central wave- 
lengths A c ~ 1.05 pm, 1.25 /im, and 1.60 /im, respectively. 
The WFC3 pointing is centered at 3 /, 32 m 38.5 s , -27°47'00." 
(J2000) and covers approximately a single WFC3 field of 
view. The total area surveyed is 4.7 arcmin 2 with maximum 
depth, taking into account dithered offsets between images. 

We processed the WFC3 data (described in more detail in 
Koekemoer et al. 2010, in prep), first reducing the data from 
the archive using the STSDAS PyRAF task calwf c3, which 
applies steps to remove the dark current, and corrects for the 
flat field and detector linearity. We used on-orbit calibration 
files for the dark current and flat-field rather than the pre- 
launch calibration files (c.f., Bouwens et al. 2009c; Oesch 
et al. 2009). Several of the F105W images suffered from 
persistence effects from a previously-observed program. We 
reduced the WFC3 images with and without these images, 
and found that it affected the F105 photometry by as much as 
several tenths of a magnitude. We therefore excluded these 
data to prevent the persistence-laden images from affecting 
our photometry. We combined the images using MultiDriz- 
zle (Koekemoer et al. 2002), aligning each image by their 
fractional-pixel dithers, and weighting by the inverse variance 
of each pixel. Each WFC3 image was aligned to the original 
ACS HUDF images to determine the unique astrometric and 
geometric transformation. Other reduction steps are similar 
to those described in McLure et al. (2009). The final mosaic 
has 0"03 pixels in each band. We measure a full-width at 
half-maximum (FWHM) of 0."176 in H\(,q for point sources 
identified in the field. 

We use the latest zeropoint information available from the 
STScI/WFC3 instrument science team 1 , which were 26.27, 

1 http://www.stsci.edu/hst/wfc3/phot_zp_lbn 
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26.25, and 25.96 AB mag for the F105W, F125W, and F160W 
images, respectively. These are identical to those employed 
by McLure et al. (2009). However, they are offset from 
those used in Bunker et al. by 0.10-0.15 mag. This offset 
is typically smaller than the statistical uncertainty in object 
^105 —^125 and .A 25 -#160 colors, however it is a systematic off- 
set in colors that can effect the interpretation of object colors, 
photometric-redshift estimation, and stellar-population prop- 
erties (see § 3 below). We measure the limiting flux depth 
of each image by measuring the distribution of fluxes in 10 4 
randomly placed 0."4-diameter apertures. We derive 5a lim- 
iting magnitudes of T105 = 28.98 mag, = 29.17 mag, and 
H\(/q = 29.21 mag. These are consistent with respect to the 
findings of other groups. 

We convolved the ACS data with a kernel to match the PSF 
of the WFC3 data. Our tests indicated that this PSF convolu- 
tion matches point-source photometry between the ACS and 
WFC3 bands to better than 5% for apertures larger that 0"5- 
diameter. Our same tests indicated that we recover the same 
fraction light to better than 5% for unconvolved versions of 
the three WFC3 images (F105W, F125W, F160W) in aper- 
tures larger than 0"5-diameter. Because the photometry mea- 
sured in the native angular resolution of these bands is already 
well matched we make no adjustment to the FWHM of the 
WFC3 images. We use the PSF-matched versions of the ACS 
data for our photometric catalogs. 

2.2. Spitzer IRAC Data 

The HUDF was observed with Spitzer at 3.6, 4.5, 5.8, and 
8.0 /Mm with the Infrared Array Camera (IRAC, Fazio et al. 
(2004)) as part of the Great Observatories Origins Deep Sur- 
vey (GOODS; Dickinson et al., in prep). The GOODS IRAC 
data are substantially deeper at 3.6 and 4.5 fim than at longer 
wavelengths, and we use only these data here. We measure 
lcr limiting magnitudes in 3.0" diameter apertures of true = 
26.94 and 26.88 mag in the 3.6 and 4.5 fim bands, respec- 
tively, using the same method as with the WFC3 data. 

Given the faint nature of the sources of interest in our sam- 
ple, they will be near the limit of the IRAC data if detected. 
Owing to the fact that the IRAC point-response function is 
substantially larger than that of HS7YWFC3 (« 0."2 FWHM 
versus 1"6 FWHM), much of the IRAC image is "contami- 
nated" by other (foreground) galaxies. Therefore, when we 
measure IRAC flux densities and upper limits in our data, we 
used a version of the IRAC image in which all sources de- 
tected in the GOODS vl.9 ACS F850LP image were removed 
using TFIT (Laidler et al. 2007). 

2.3. Object Detection and Photometry 

We used the Source Extraction (SExtractor) software pack- 
age v2.5.0 (Bertin & Arnouts 1996) in two-image mode to 
construct an object catalog and perform photometry. We con- 
structed a detection image using a sum of the F125W and 
F160W images weighted by the image inverse variance fil- 
tered by a 3-pixel Gaussian kernel. Objects were detected at 
> 2cr significance with a minimum area of 12 pixels. Objects 
were then photometered on each of the ACS F435W, F606W, 
F775W, F850LP, and WFC3 F105W, F125W, and F160W im- 
ages. We measure object photometric colors using elliptical 
apertures that scale with the size of each source defined fol- 
lowing Kron (1980), with a factor (KRON_FACT) of 1.2 and a 
minimum size (MIN_RADIUS) of 1.7. Our tests have shown 
that for the faint galaxies in our sample {H\^q > 28.5 mag) 



these Kron-like apertures yield fluxes comparable to those 
measured in either isophotal apertures or 0."4-diameter aper- 
tures. For brighter galaxies in our sample (H < 28.0 mag), 
these Kron-like apertures contain substantially more flux, in 
some cases by as much as 1 magnitude. Furthermore, we 
find that the signal-to-noise ratios (S/N) in these Kron-like 
apertures are comparable to those measured in isophotal aper- 
tures, but because they are based on the light profile of each 
galaxy they produce higher S/N than that measured in circu- 
lar apertures by 10-30% for all galaxy magnitudes. Therefore 
we adopt the Kron-like apertures to measure galaxy colors as 
they maximize the S/N while containing an optimal light frac- 
tion. To measure galaxy total fluxes, we apply an aperture cor- 
rection to each band constructed from the weighted-summed 
/i25 + fli6o detection image. The aperture corrections are de- 
fined on a source-by-source basis as the difference between 
the light in the Kron-like aperture defined above to a larger 
Kron-like aperture, with a Kron factor of 2.5 and a minimum 
radius of 3.5. 

2.4. Photometric Redshifts 

We derive redshift estimates for each galaxy in our WFC3- 
based catalog using the full ACS B435, V606, '775, zsso an d 
WFC3 T105, /125, #160 photometry and detections or upper 
limits from IRAC [3.6], [4.5]. For the ACS and WFC3 pho- 
tometry, we use flux densities measured for each galaxy in the 
Kron-like elliptical aperture defined in § 2.3. We scale these 
to total magnitudes using the aperture corrections defined in 
§ 2.3. For the IRAC data, we use the measured flux densities 
from the GOODS IRAC catalog (Dickinson et al., in prep) 
for those objects with detections. For the remaining sources, 
we use upper limits measured from the image "cleaned" of 
sources detected in the GOODS z%5Q image, as described in 
§ 2.2. We correct the IRAC fluxes for light falling outside 
the circular apertures by applying corrections of 0.36 mag and 
0.3 1 mag (derived for point sources) to the 3.6 /im and 4.5 fim 
data, respectively. The upper limits are handled identically to 
detections, i.e. the measured flux in the aperture is used in the 
fitting, and the 1 a error is the noise. The only exception is 
when the measured flux is formally negative - in these cases, 
the flux is set to zero, and the 1 a error is still used as the 
noise. 

We use the photometric-redshift code EAZY (Brammer 
et al. 2008) to estimate photometric redshifts and derive red- 
shift probability distribution functions, P(z), for the sources in 
our WFC3 catalog. The measured photometry for each galaxy 
was fit with non-negative linear combinations of a suite of 
optimized templates provided with EAZY, based on the PE- 
GASE stellar-synthesis models (Fioc & Rocca-Volmerange 
1997). We also include the composite rest-frame UV SED 
for z ~ 3 LBGs from Shapley et al. (2003) as this may best 
represent the galaxies of interest in our sample and it con- 
tains Lyman a emission. We included the effects of absorp- 
tion from intervening H I in the intergalactic medium (IGM) 
along the line of sight following the prescription of Madau 
(1995), though we note the Lya forest is effectively opaque at 
these redshifts. 

We fit photometric redshifts to all objects in our (/125 + 
i/i6o)-selected catalog, including the 4 ACS bands, 3 WFC3 
bands, and the 2 bands from the IRAC data. For the purposes 
of photometric-redshift analysis, we included those galaxies 
in our catalog with /125 < 29.55 mag and H\gq < 29.60 mag, 
which are equivalent to the 3.5cr-detection limits for sources 
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FIG. 2. — 5"cutouts of our 31 candidate z > 6.3 galaxies, from left-to-right: (smoothed) ACS i' 115 and Zg 5Q , and WFC3 Y105, J125 and ffi6o- The objects 
are shown in order of increasing photometric redshift. The circles are centered on the objects and have a radius of 1". The last column shows the probability 
distribution function for the photometric redshifts, and prints the integrated probability over 6 < z < 1 1 CPs). 



in 0."4-diameter apertures. For each galaxy we calculate the 
photometric redshift with the lowest x 2 , as well as the full 
photometric redshift probability distribution function, P(z), 
defined as P(z) <x exp(~x 2 /2), normalized such that the in- 
tegral of P(z) over all redshifts is equal to unity. 



2.5. Selection of Galaxies at 6.3 < z < 8.6 

We used the full redshift probability distribution function, 
P(z), defined above to select those targets with a high like- 
lihood of having redshift z > 6. We define the "integrated 
probability", Vd, as the integral of the redshift probability dis- 



tribution function from 6 < z < 11, that is 



V 6 = J P(z)dz. 



(1) 



For our sample, we selected all objects satisfying Vb > 0.6, 
thus each of these galaxies has >60% of its probability distri- 
bution function above redshift z = 6. This method excludes no 
information on the galaxies. We make no other requirements 
on any of the galaxy colors nor on the limits (non-detections) 
in any of the ACS bands. This is in contrast to other typi- 
cal UV-"dropout" techniques targeting the redshifted Lyman 
break (see Bouwens et al. 2009c; Oesch et al. 2009; Bunker 
et al. 2009; Yan et al. 2009, and discussion below). Instead, 




our photometric redshifts use the full set of photometry from 
all bandpasses including all of the photometric uncertainties 
to construct the likelihood that a galaxy lies at a given red- 
shift, although this method primarily keys off of the Lyman 
break feature. 

Our sample of galaxies satisfying /125 < 29.55 mag, H^q < 
29.60 mag and > 0.6 includes 45 galaxies. Upon visual 
inspection of these sources, we exclude 3 objects as obvious 
spurious sources (stellar diffraction spikes, oversplit areas of 
bright foreground galaxies). We also exclude one additional 
source with /125 = 26.2 mag and Hk,o = 26.1 mag. This object 
has been noted by others (Oesch et al. 2009; McLure et al. 
2009; Yan et al. 2009) as it is bright, yet undetected in ex- 
isting NICMOS F110W and F160W imaging of the HUDF, 
which implies it is a transient source. Furthermore, this ob- 



ject is unresolved at the WFC3 resolution, with a measured 
FWHM=0"18 (compared to Of 176 PSF FWHM for WFC3, 
see also discussion below). Our final sample of high-redshift 
objects includes 41 galaxies. 

Figure 1 shows the number distribution of the best-fit 
photometric-redshifts, z p hot, of the galaxies in this sample, 
which have P(z) distributions that are sharply peaked above 
Z = 6. From this parent sample, we construct two redshift 
samples for study here denoted as the "Sample at z ~ 7" and 
"Sample at z ~ 8" for convenience. These samples are defined 
as, 



Sample at z ~ 7 : 6.3 < z p hot < 7.5 
Sample atz ~ 8 : 7.5 < z p hot < 8.6. 
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FIG. 2. — Continued. 
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The z ~ 7 sample includes 23 galaxies 2 with a mean red- 
shift (z) = 6.73 and the z ~ 8 sample includes 8 galaxies with 
(z) = 8.00. For our analysis, we will focus primarily on these 
31 galaxies with 6.3 < z p /, or < 8.6. Figure 1 shows also the 
joint redshift probability distribution functions for each sam- 
ple, defined as Pj j n t(z) = Sf,(z), where the Pj{z) are the red- 
shift probability distribution functions for the ith galaxy in 
each sample. We have tested that the z ~ 7 and z ~ 8 subsam- 
ples are not subject to our arbitrary definition of the integrated 
probability distribution function, Ve- Lowering the integra- 
tion limits on V from zi OW ei- = 6 to zi 0W er = 5 has no change on 

2 One of the galaxies in the z = 1 sample (ID 1441) was split into two 
separate objects by SExtractor. Based on the similarity in the object colors 
and photometric redshift probability distribution functions, we have merged 
these objects by hand in our sample. 



the objects in these samples. Furthermore, the majority of the 
objects in these samples have high . Of the 3 1 objects, 27 
(23) have V b > 0.8 (0.9). We chose the limit at V 6 > 0.6 to 
maximize our sample size while minimizing the contamina- 
tion rate, which we found becomes higher at V§ < 0.6 based 
on sources not satisfying the LBG photometric criteria (see 
§3.1). Out of our 31 galaxies, 4 have z < 5 contained within 
the 68% confidence range of their photometric redshift prob- 
ability distributions, yet 3/4 satisfy the zsso- or Y\ 05 -dropout 
criteria (Oesch et al. 2009; Bouwens et al. 2009c). To keep 
the interpretation physical, we applied a Bayesian prior re- 
stricting the redshift to > 6.5 or 7.5, respectively, for these 
three objects. The remaining object, ID 1818, retains its large 
redshift uncertainty. 
Table 1 provides the astrometric coordinates, photometry, 
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and photometric-redshift information of the 23 objects in 
the z ~ 7 sample and the 8 objects in the z ~ 8 sample. 
Figure 2 shows postage-stamp images in the ACS F775W, 
F850LP and WFC3 F105W, F125W, and F160W bands, and 
the P(z) distribution, for each of the galaxies in the z ~ 7 and 
z ~ 8 samples. Generically, we will designate the galaxies 
in these samples as LBGs. Although they were not selected 
strictly based on the strength of their Lyman breaks, this is 
the key feature driving their photometric redshifts, and as a 
result there is substantial overlap with LBG-selected samples 
(see § 3.3). Each galaxy in our sample is formally unde- 
tected in the B435, V606, and 2775 bands (with S/N<2.0 in each 
band, with the exception of ID 1566, which has S/N, = 2.7 3 ), 
and they each exhibit large color indexes between either the 
Z850 _ i'los, or F105 - 7i25 bands, indicative of the Lyman-break 
redshifted into these bandpasses. 

2.6. IRAC Photometry of Candidate High-Redshift Galaxies 

3 Due to image convolution, some flux from an adjacent galaxy bleeds into 
this object's (775— band aperture. 



In order to measure the IRAC flux at the position of each of 
our candidates, we constructed a simulated detection image 
using the IRAF 4 task mkob ject s in the art data package, 
populated with sources with Gaussian profiles inserted at the 
positions of the high-redshift sources in our sample. We then 
performed photometry in 3.0" diameter circular apertures on 
the source-subtracted IRAC images using SExtractor in two- 
image mode. 

While Labbe et al. (2009a) detect z ~ 7 galaxies with IRAC, 
their objects were selected from the shallower WFC3 Early 
Release Science (ERS; Windhorst et al. 2010) data, and are 
thus intrinsically brighter. In contrast, only one object in 
our high-redshift sample (ID 2013) is detected in the IRAC 
3.6 /.tin data at > 3 a significance (3.29 a). Upon visual in- 
spection of this source, we found it to be strongly blended 
with an adjacent galaxy, with large residuals in the TFIT- 

4 The Image Reduction and Analysis Facility (IRAF) is distributed by the 
National Optical Astronomy Observatory (NOAO), which is operated by the 
Association of Universities for Research in Astronomy, Inc. (AURA) under 
cooperative agreement with the National Science Foundation. 
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TABLE 1 
Candidate Galaxies at z > 6.3 



ID 


RA 


Dec 


Zphot 


v 6 


j 

z 850 


^105 


Jl25 


#160 


*?150O 


SFR,/v Flags 




(J2000) 


(J2000) 








(hub) 




(m A s) 


(mag) 


(Ms yr 1 ) 



803 


53. 


18242 


-27 


.77603 


2894 


53, 


15322 


-27 


.79822 


769 


53, 


15610 


-27 


.77578 


1289 


53. 


15601 


-27 


.78091 


649 


53. 


17579 


-27 


.77440 


1445 


53, 


16491 


-27 


.78235 


2032 


53. 


15158 


-27 


.78784 


1818 


53. 


17601 


-27 


.78576 


328 


53, 


17427 


-27 


.76979 


1072 


53. 


16903 


-27 


.77876 


1566 


53. 


14554 


-27 


.78372 


2013 


53. 


17104 


-27 


.78767 


1110 


53, 


18642 


-27 


.77913 


1441 


53. 


17735 


-27 


.78240 


567 


53. 


17677 


-27 


.77342 


2432 


53, 


17734 


-27 


.79206 


515 


53. 


16553 


-27 


.77259 


335 


53, 


15984 


-27 


.76997 


1768 


53, 


16169 


-27 


.78532 


2056 


53. 


16481 


-27 


.78819 


669 


53, 


17973 


-27 


.77457 


1106 


53, 


18626 


-27 


.77897 


3053 


53. 


15506 


-27 


.80171 


819 


53. 


17866 


-27 


.77625 


653 


53, 


17952 


-27 


.77436 


2055 


53. 


16466 


-27 


.78815 


200 


53, 


15749 


-27 


.76670 


213 


53, 


15681 


-27 


.76709 


3022 


53. 


16806 


-27 


.80073 


640 


53. 


15334 


-27 


.77457 


125 


53, 


15890 


-27 


.76500 



6 34+0.17 
-0 22 
6 14 rt:i! 
6 3gai.n 

6.J8_ 016 
6 40+4). 18 

6 40^^ 

640 +4llS 

6.43+0.31 
6 47+<«J 
6 47+0-25 
6 56+* 
6.70 rt33 , 
6 72+^ 

b - SJ -0.13 
6 - 86 l 

6.887$ 

6 93+* 37 
y 22+0-18 

7 25 +8:M 

'■ ZJ -0.23 
7.25+ 032 

7 o9+8:55 

-0.20 
7 40^- 3ii 

7.55^- 3 i 



7.76 
7.81 



-0.33 

7 91+4)19 

8 05+^ ^ 

8.24+0- 3 ' 
8.61*| 



0.93 


28.96 ± 0.29 


28.04 ±0.13 


28.15 ±0.12 


28.17 ±0.11 


-18.71 


13.7 


3 


>0.99 


28.61 ±0.21 


27.76 ±0.10 


27.85 ± 0.09 


27.96 ± 0.09 


-19.00 


17.8 


2.3 


0.96 


27.38 ± 0.09 


26.68 ± 0.05 


26.64 ± 0.04 


26.47 ± 0.04 


-20.12 


49.9 


3 


0.99 


28.89 ± 0.29 


27.81 ±0.11 


27.98 ±0.11 


28.05 ±0.11 


-18.94 


17.0 


2,3,5 


>0.99 


28.53 ± 0.24 


27.69 ±0.11 


27.53 ± 0.08 


27.73 ± 0.09 


-19.16 


20.8 


3 


0.96 


> 29.24 


28.23 ±0.15 


28.23 ±0.12 


28.68 ±0.18 


-18.57 


12.1 


2,3,5 


0.99 


28.97 ± 0.36 


27.74 ±0.12 


28.00 ±0.13 


28.05 ±0.13 


-18.99 


17.8 


2,3,4,5 


0.76 


> 28.28 


27.66 ± 0.21 


28.07 ± 0.25 


28.02 ± 0.23 


-19.03 


18.5 


- 


0.89 


> 29.46 


28.85 ± 0.21 


28.89 ±0.18 


29.05 ± 0.21 


-17.95 


6.9 


5 


0.98 


> 29.11 


28.15 ±0.15 


28.29 ±0.15 


28.47 ±0.16 


-18.61 


12.7 


2,3,4,5 


>0.99 


28.53 ± 0.21 


27.27 ± 0.07 


27.17 ±0.05 


27.31 ±0.06 


-19.61 


32.2 


3 


>0.99 


> 29.18 


28.42 ±0.18 


28.05 ±0.11 


28.49 ±0.16 


-18.66 


13.7 


2,3,5 


0.80 


> 28.39 


28.20 ±0.31 


27.94 ± 0.20 


27.93 ±0.19 


-18.82 


15.9 


3,5 


>0.99 


27.77 ±0.16 


26.13 ±0.04 


25.95 ± 0.03 


25.95 ± 0.02 


-20.89 


108.1 


2,3,4,5 


0.88 


> 28.63 


28.43 ± 0.30 


28.41 ± 0.25 


28.57 ± 0.28 


-18.48 


11.8 


5 


>0.99 


> 28.93 


27.38 ± 0.09 


27.17 ±0.06 


27.25 ± 0.06 


-19.66 


35.2 


2,3,4,5 


0.99 


> 29.15 


28.56 ± 0.21 


28.46 ±0.16 


28.70 ±0.19 


-18.42 


11.2 


2,3,4,5 


0.99 


> 29.55 


28.90 ± 0.20 


28.47 ±0.11 


28.65 ±0.13 


-18.34 


10.5 


4,5 


>0.99 


> 29.16 


27.50 ± 0.08 


27.07 ± 0.05 


27.01 ± 0.04 


-19.89 


45.2 


2,3,4,5 


>0.99 


> 28.71 


27.72 ±0.14 


27.15 ±0.07 


27.05 ± 0.06 


-19.82 


42.5 


2,3,4,5 


>0.99 


> 28.76 


28.03 ±0.19 


27.58 ±0.10 


27.50 ± 0.09 


-19.39 


28.8 


2,3,4,5 


>0.99 


> 29.01 


27.87 ±0.13 


27.27 ± 0.06 


27.45 ± 0.07 


-19.73 


39.7 


2,3,4,5 


>0.99 


> 29.04 


28.45 ± 0.21 


27.77 ± 0.09 


27.96 ±0.11 


-19.24 


25.5 


2,3,4,5 


>0.99 


> 29.10 


28.82 ± 0.28 


28.09 ±0.12 


28.11 ±0.12 


-18.96 


20.1 


1,3,4,5 


0.95 


> 29.06 


> 29.04 


28.31 ±0.15 


28.92 ± 0.26 


-18.68 


15.8 


1,3,5 


0.96 


> 28.64 


28.41 ± 0.30 


27.58 ±0.11 


27.21 ± 0.08 


-19.63 


38.3 


2,3,4,5 


0.89 


> 29.06 


> 29.05 


28.02 ±0.12 


28.29 ±0.14 


-19.03 


22.2 


1,3,4,5 


0.86 


> 29.32 


> 29.32 


28.70 ±0.17 


28.77 ±0.18 


-18.43 


12.9 


1,3,4,5 


0.62 


> 29.36 


> 29.35 


29.15 ±0.26 


28.97 ± 0.21 


-18.07 


9.3 




0.62 


> 27.67 


> 27.62 


27.23 ± 0.20 


26.96 ±0.15 


-20.10 


61.6 




0.95 


> 29.48 


> 29.47 


28.61 ±0.14 


28.31 ±0.11 


-18.85 


20.4 


1,3,4,5 



NOTE. — The objects are listed in order of increasing photometric redshift, and the magnitude upper limits shown are 3 a. The 
photometric redshift errors are 1 a. The column Ve is the integrated probability from 6 < z < 11, where only objects with Vb > 0.6 
were included in our sample. M1500 denotes the absolute magnitude at rest-frame 1500 A, which was computed for each object by 
converting the WFC3 fluxes to rest-frame (with z = Zphot), and interpolating to 1500 A. The SFRs were computed using M1500 along with 
the Kennicutt (1998) UV luminosity density - SFR relation. The flag column denotes which objects were previously discovered by other 
studies, with the numbers representing: 1) Bouwens et al. (2009c); 2) Oesch et al. (2009); 3) McLure et al. (2009); 4) Bunker et al. (2009) 
& 5) Yan et al. (2009). Only ID 2056 was detected at > 3 a in either of the IRAC bands, with m [4 .5] = 25.53 ±0.31. The remaining 
objects have 3 a upper limits of > 25.75 in [3.6] and > 25.69 in [4.5]. 



subtracted image. We therefore discard this flux measurement 
and replace its IRAC flux with zero in our subsequent anal- 
ysis, using the measured 3.6 fim magnitude limit as the la- 
error. 

Four objects in our high-redshift sample are detected for- 
mally in the IRAC 4.5 /im image, with 2.6-3.5 a significance. 
None of these sources were detected in the IRAC 3.6 /.tm im- 
age. Two of the objects with 4.5 /im-detections (ID 335 & 
1 8 1 8) are blended with light from galaxies that are poorly sub- 
tracted by TFIT, and we set their 4.5 /im flux to zero, similar 
to how we handled the 3.6 fim source discussed in the previ- 
ous paragraph. 

The other two objects, source IDs 2055 and 2056, are de- 
tected at 4.5 fim with 3. lcr and 3.5cr significance, respectively. 
Based on the photometric redshift estimates of these objects 
(see § 2.4) it is possible that the IRAC 4.5 is contami- 
nated by emission from [O III] AA4959,5007, providing a pos- 
sible explanation for the 4.5 fim detections and lack of flux at 
3.6 fim. However, these objects lie within ~ 1" from each 



other, with overlapping isophotes in the WFC3 images, and 
they are blended at the resolution of IRAC. 5 . Furthermore, 
we find that each of these galaxies have photometric red- 
shift probability distribution functions consistent within their 
68% confidence ranges (see § 2.4), with z p h t = 7.8 l^P> and 
7.25^Q93, respectively. Given their close proximity and simi- 
lar redshifts (within the uncertainties) , these objects may be 
physically associated. However, because we are unable to de- 
blend the IRAC data, we assign all the 4.5 fim flux to object 
ID 2056 because this object is brighter by 0.2 mag in i/]6o and 
brighter by 0.4 mag in /125. Even in the case that they are are 
unassociated physically, ID 2056 would likely dominate the 
emission. 

3. COLOR-SELECTION OF GALAXIES AT 6.3 < z < 8.6 

5 The photometry at the positions of these two objects was performed sep- 
arately, thus the detection significance of the combined IRAC fluxes from 
both is still ~ 3 <j. 
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FIG. 3. — Color-color (Fios - Jns versus z'g 50 ■ 
photometric redshift selected objects compare to the selection criteria from 
various studies. Small gray squares denote objects with z p h t < 6.0. Filled 
gray circles denote objects with 6.0 < Zphot < 6.3. Filled cyan circles denote 
objects with 6.3 < z p hot < 7.5. The arrows denote 1 a limits. The solid lines 
show the selection criteria from Oesch et al. (2009), while the dashed line 
shows the selection criterion from Yan et al. (2009), both for LBGs at z ~ 7. 
Due to our photometric redshift analysis our sample includes galaxies outside 
of these windows. Brown stars show the expected colors of galactic brown 
dwarfs. Objects inside the LBG selection window which are not selected as 
high-redshift galaxies have significant optical flux, resulting in low-redshift 
photometric redshift fits. 



3.1. The Colors of Galaxies at 6.3 < z < 8.6 

Figure 3 shows the T105-./125 versus zsso-^ios color-color 
plot for galaxies in the HUDF, which is useful to isolate z ~ 7 
LBGs. Galaxies with z > 6.3 reside in the upper-left portion 
of the color-color space. These objects have red z^sq - Y W5 
colors, indicative of the Lya break redshifted to A > 8900 A. 

We also show the color-selection criteria of Oesch et al. 
(2009) and Yan et al. (2009) used to select z ~ 7 galaxies 
in Figure 3, in order to contrast our selection method us- 
ing the full photometric redshift probability distribution func- 
tion with pure color-selected samples. While these classical 
color-selection criteria identify the majority of objects in our 
sample, there are distinct differences. Our sample selected 
using the full photometric redshift P(z) includes galaxies ex- 
cluded by these color-selection criteria. Parenthetically, we 
find similar results if we use the photometric-redshift-derived 
samples from McLure et al. (2009). Because both the crite- 
ria of Oesch et al. and Yan et al. require Z850 _ Tio5 > 0.8 mag 
for object selection, they miss some objects with colors just 
below this threshold (owing to a combination of photometric 
uncertainties and color variations). This is also expected as 
galaxies with intrinsically bluer F105— J125 colors enter these 
color-selection windows at lower redshifts than galaxies with 
redder T105-/125 colors. This is illustrated by the fact that 
there is a source with T105-/125 = -0.3 mag and z p h t < 6.3 
within the Oesch et al. color-selection "wedge", whereas sev- 
eral galaxies with z p h t > 6.3 with redder T105-/125 colors lie 
outside the color-selection windows. 

Brown-dwarf stars are expected to have red Z&50 - T105 col- 
ors, mimicking the colors of high-redshift galaxies. Fig- 
ure 3 shows the synthesized colors of brown-dwarf stars (us- 




>- 1 
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1.0 



0.0 

J 125 — H 160 

FIG. 4. — 7j25 - Hi(,o versus Fios - /125 color-color plot for galaxies in 
our WFC3 sample. Symbols are the same as in Figure 3, except here we 
also show galaxies with 7.5 < Zphot < 8.6 as red filled circles. The solid lines 
show the selection criteria from Bouwens et al. (2009c), while the dashed line 
shows the selection criterion from Yan et al. (2009), both for z ~ 8 LBGs. 

ing spectra of brown dwarfs from S. Leggett 6 , with the data 
from Golimowski et al. (2004), Knapp et al. (2004) and Chiu 
et al. (2006)). All the brown dwarfs in this sample would sat- 
isfy Z850 _ Tio5 > 0.8 mag. The Oesch et al. color-selection 
"wedge" would exclude most of these objects. However, all 
of these objects would be present in a sample using the cri- 
terion of Yan et al., although the authors argue their sample 
includes no unresolved sources. 

Figure 4 shows the /125 -#i60 versus T105 -^125 color-color 
plot for galaxies in the HUDF, which is useful to isolate z ~ 8 
galaxies. Those galaxies in our z ~ 8 sample are indicated by 
red-filled circles, while galaxies in the lower-redshift z ~ 7 
sample are indicated as cyan-filled circles. Similar to the z ~ 
7 sample in Figure 3, the z ~ 8 sample resides largely in the 
upper-left portion of plot, and they have red T105 - /125 colors, 
indicative of the Lya break redshifted to A > 1 1000 A. 

Figure 4 also shows the color-selection criteria of Bouwens 
et al. (2009c) and Yan et al. (2009) to select z ~ 8 galaxies. Of 
the eight galaxies in our z ~ 8 sample, one with T105 — /125 = 
0.7 mag (ID 819) would be missed by both color-selection 
criteria, and one additional galaxy (ID 2055) with T105 — /125 = 
0.8 mag and /125 - H^q = 0.4 mag lies outside the Bouwens 
et al. "wedge". 

As with Figure 3 brown-dwarf stars provide contamination 
to z ~ 8 samples. Figure 4 shows that the expected /125 — #160 
colors of brown dwarfs span a large range. While most have 
Tio5 - -/i25 < 0.8 mag, the distribution lies near to this limit, 
implying that photometric uncertainties of brown dwarfs will 
mimic the colors of z ~ 8 galaxies. 

3.2. The Light Profiles of Galaxies at 6.3 < z < 8.6 

Owing to the fact that the expected Zgso-Tios, T105-/125, 
and 7i25 -//16O colors of brown-dwarf stars can mimic those 
of z > 6 galaxies, we examined the size of the light profiles of 
the galaxies in our sample. Figure 5 shows the full-width at 
half-maximum (FWHM) of the galaxies in our catalog mea- 

6 http://staff.gemini.edu/~sleggett/LTdata.html 
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sured from their light profiles in the WFC3 F160W image. 
The figure indicates FWHM = 0."18, the expected FWHM of 
WFC3 PSF in the F160W band. We identified all bright stars 
in the WFC3 F160W image, and these all have a FWHM w 

6 pixels (0."18 at 0."03 pix" 1 ). We subsequently examined all 
objects in our full WFC3 catalog with a FWHM < 6.6 pix- 
els and found them to have stellar morphologies, consistent 
with point sources. These sources are indicated as large filled 
pentagrams in Figure 5. The temporally transient source iden- 
tified originally in our sample of high-redshift galaxies is also 
indicated in the figure, and it is consistent with being a point 
source. In contrast, the FWHM measured for the F160W light 
profiles of our high-redshift samples are resolved, and exceed 
the FWHM of point sources down to the magnitude limit of 
our catalog. From this we conclude that the sources in our z ~ 

7 and z ~ 8 samples are resolved galaxies, with no contami- 
nation from stellar sources (including brown dwarfs). 

3.3. Comparison to High-Redshift Sample Selection from the 

Literature 

In this section we compare the galaxies in our z ~ 7 sample 
and z ~ 8 sample to others in the literature using these WFC3 
HUDF data. Table 1 provides a flag for sources in our sam- 
ple that match the astrometric coordinates of sources in these 
other studies. In the subsections below we provide qualitative 
comparisons. 

3.3.1. Comparison to Bouwens et al. & Oesch et al. 

Bouwens et al. (2009c) examined an early reduction of the 
HUDF WFC3 data, searching for "Y\ 05 "-dropout galaxies at 
7.5 < z < 8.5 (see also Figure 4). Using these criteria with 
additional criteria to reject low-redshift interlopers, they re- 
ported five candidate LBGs, all with H160 > 28 mag. We re- 
cover all five of these z ~ 8 candidate LBGs from Bouwens 
et al. (2009c), and we find that the photometric redshifts of 
these galaxies range from 7.6 < z p hot < 8.6. The measured 
J\25-H\(,q colors of the matched objects are consistent within 
< 0.2 mag. 

Oesch et al. (2009) used the same reduction of the WFC3 
HUDF data as Bouwens et al. and searched for "z85o"-dropout 
galaxies at 6.5 < z < 7.5 (see also Figure 3). Oesch et al. iden- 
tified 16 Zg 50 -dropout galaxies. We match 14 of these galaxies, 
for which we derive a photometric redshift range from 6.3 < 
Zphot < 7.4. Of the two objects not in our catalog, one fails 
our magnitude limit for selection as it has /125 = 29.6 mag and 
#160 = 30.3 mag, and for the other we compute z p hot = 5.7. The 
■A25 _ #160 colors we measure for these galaxies are consistent 
within < 0.2 mag for 12 out of the 14 objects in common be- 
tween our catalogs. For the remaining objects, we measure 
color discrepancies as high as 0.28 mag, which we attribute to 
variations in data reduction (see § 2.1). 

3.3.2. Bunker et al. 

Bunker et al. (2009) perform color-color selections similar 
to that employed by Oesch et al. (2009) and Bouwens et al. 
(2009c), selecting 1 1 z g50 -dropout galaxies with z ~ 7, and 
7 Fios-dropout galaxies with z ~ 8. We match all 1 1 Zg 50 - 
dropouts from Bunker et al. and we measure a photometric- 
redshift range of 6.4 < z p h t < 7.4. We note, however, that 
in our catalog object zD4 of Bunker et al. (2009) is split into 
two objects in our sample, ID 2055 and 2056 (see §2.6). We 
recover 4 of the 7 T^-dropouts from Bunker et al., measuring 
a range of 7.5 < z p hot < 8.6. We excluded 2 of the missing 
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FIG. 5. — Because Galactic brown dwarfs can have similar colors to high 
redshift galaxies, we compare the measured FWHM of our sources to stars 
in the image. Unresolved sources at the WFC3 resolution have FWHM < 6 
pixels (Of 18). We visually inspected unresolved objects in the H^q image, 
and all appear very compact and consistent with being stellar. We note that 
one of these objects is the transient first detected by Oesch et al. (2009). All 
of our candidate high-redshift galaxies are resolved, with FWHM > 0"2. We 
conclude our sample suffers no contamination by galactic brown dwarfs. 

3 sources selected by Bunker et al. because they are fainter 
than our 3.5er magnitude limit in H\sq. The remaining source 
(Bunker et al. ID YD5) has a photometric redshift of z p hot = 
2.0^Q g, and was therefore excluded from our sample. 

3.3.3. McLure et al. 

The selection of McLure et al. (2009) is most similar to that 
used here. McLure et al. used an independent reduction of the 
WFC3 HUDF data, and selected all WFC3 sources with no 
detected counterpart in the ACS F775 W images. We match 25 
of the 35 objects identified by McLure et al. (2009) at z > 6.3. 
Of the 10 objects not in our catalog, 5 have 5.7 < z p h t < 
6.2, one has z p hot = 1 .SIq'I, and 3 are fainter than our catalog 
magnitude limits. The remaining object has z p h t = 6.6^54, 
with the large redshift uncertainty giving V(, < 0.6. 

For the 25 objects in common between our catalogs, we find 
that our photometric redshifts are consistent within Az ss 0.2 
for 17 objects. Although our /125 - //16O colors are consis- 
tent with that measured by McLure et al. (< 0.3 mag for all 
sources), small differences due to photometric aperture vari- 
ations affect the redshift interpretation. The catalog used by 
McLure et al. used 0"6-diameter circular apertures to mea- 
sure colors, in contrast to the elliptical apertures used to con- 
struct our catalog. 

3.3.4. Yanetal. 

Yan et al. (2009) used an independent reduction of the 
WFC3 HUDF data to select z 850 -dropout (z ~ 7), Ti 05 - 
dropout (z ~ 8), and /125-dropout (z ~ 9) galaxies. They iden- 
tify 20 z 850 -dropouts, 15 Tios-dropouts, and 20 /125 -dropouts. 
We match 16 of the 20 z 850 -dropouts in our catalog, for which 
we measure a photometric redshift range of 6.4 < z p h t < 
7.4. We match 7 of the 15 Y\ 05 -dropouts, with 6/7 residing 
in our z ~ 8 sample, and one additional having z p h t = 6.8 (in- 
cluded in our z ~ 7 sample). Of the 12 galaxies of Yan et 
al. unmatched in our samples, 8 are fainter than the magni- 
tude limit required for our selection, including 1 galaxy that 
is split into two components in our catalog (both of which are 
then fainter than our magnitude limit). For three others we 
derive photometric redshifts z p hot < 6.3. We excluded the re- 
maining galaxy because it has overlapping isophotes with a 
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bright nearby galaxy, which appears to contaminate its col- 
ors. The additional galaxies in the sample of Yan et al. are 
fainter and less secure as Yan et al. require only a 3cr detec- 
tion in either the F125W or F160W-bands. Additionally, Yan 
et al. did not restrict the /125 -#i60 color of their I 7 ! 05 -dropout 
sample, which possibly suffers a higher contamination from 
galactic brown dwarfs (although Yan et al. exclude unresolved 
sources) and lower redshift galaxies with older stellar popu- 
lations (see also Figure 4). As detections were required in 
both the F125W and F160W -bands for our photometric red- 
shift analysis, we have no sample to compare with the single- 
band-detected 7i25 -dropouts of Yan et al. 

4. REST-FRAME UV PROPERTIES OF GALAXIES AT 6.3 < z < 8.6 

4.1. Rest-Frame UV Colors 

Figure 6 shows the color-magnitude diagram of H^q ver- 
sus J\25-H\a) for our z ~ 7 sample. At these redshifts the 
Ji25~Hia) color probes rest-frame 1500-1900 A. This region 
is highly sensitive to massive star formation and the amount of 
dust attenuation. Even though the photometric uncertainties 
on individual objects are large owing to their faint observed 
fluxes, the z ~ 7 galaxies have blue /125 -#i60 colors. In few 
cases, these blue colors are consistent with the expected rest- 
frame UV colors of local starbursts (Kinney et al. 1996), but 
the majority of the galaxies are bluer than these local star- 
bursts. This is not a consequence of our selection, as we made 
no restriction on rest-frame UV color. To verify this, in Fig- 
ure 6 we shade the region excluded by our magnitude limits 
{J\25 < 29.55 and H^q < 29.60). This region would exclude 
galaxies only with J^s-Hiw redder than local starbursts at 
ffieo < 29. 

To quantify the photometric scatter in our data, we ran sim- 
ulations whereby we generated ^5 x 10 4 mock galaxies and 
placed them in the WFC3 data. These mock galaxies were 
generated with a uniform distribution of luminosities and red- 
shifts over 6 < z < 9. We used models from the 2007 update to 
Bruzual & Chariot (2003; hereafter CB07) for their intrinsic 
SEDs, with a constant star formation history, Salpeter initial 
mass function (IMF), age of ~ 10 8 yr and Z = 0.005 Z Q . They 
were attenuated by dust via the extinction curve of Calzetti 
et al. (2000) with -0.3 < E(B-V) < 0.5 (where the negative 
color excesses simulate extremely blue galaxies) as well as 
by the IGM (Madau 1995). The galaxy shapes were chosen 
to be exponential disks, with a log-normal radius distribution 
over the range of 0."03 < r < 1.0", with the mean recovered 
galaxy radius of 0. 14". These mock galaxies were then placed 
in the WFC3 data with the appropriate amount of noise using 
art data and we extracted their flux in much the same way 
as with our observations. Upon investigation of these cata- 
logs, we found that while the offset between the input color 
and measured color is small (< 0.05 mag from Hk,q = 25-29 
mag), the photometric scatter can be large, from a = 0.05 mag 
at //j6o = 25, to ~ 0.2 mag at H\(,q = 29. This scatter is repre- 
sentative of the true photometric uncertainty on the colors of 
our objects. We note that the scatter of ~ 0.2 mag at m# = 29 
is similar to what we obtained from our photometry in §2. 

Given the large photometric scatter on individual galax- 
ies, we compute the mean 7125 - //160 color as a function 
of Hk,o brightness, in subsamples of Z/160 > 28.3 mag and 
//160 < 28.3 mag (where = 28.3 is equivalent to Myy = 
-18.75 assuming z = 7 and the mean /125 — #160 color of the 
z ~ 7 sample). We compute the mean and uncertainty using 
Monte Carlo bootstrap simulations of 10 7 realizations of the 



data. In each simulation, we perturb randomly the colors of an 
object by its associated photometric scatter, and then compute 
the median /125 — i/160 of the simulated objects. Because these 
sample sizes are small (i.e. the H\(& > 28.3 subsample at z ~ 7 
consists of 7 galaxies) we include the effects of Poisson noise. 
We do this by including N + Axa/A 7 galaxies in the median 
for a given simulation, where N is the number of galaxies in 
a given sample, A is a random number drawn from a normal 
distribution with a mean of zero and a standard deviation of 
one, and the galaxies chosen are randomly self-sampled. 

The mean and standard deviation of this distribution is thus 
an estimate of the true mean and uncertainty on the /125 -#160 
distribution of the real sample. The histograms in the right 
panel of Figure 6 shows the results of our simulations for the 
z ~ 7 sample, and Table 2 gives the measured mean and 68% 
confidence range. 

Figure 7 shows the H\at versus /125 -#160 color-magnitude 
diagram for the z ~ 8 sample. We performed a similar calcu- 
lation of the mean and 68% confidence range for this sample, 
using the same Monte Carlo bootstrap method as for the z ~ 

7 sample discussed above. At z ~ 8 the Lya emission line re- 
sides in the F125W-band, and we indicate the its effect on the 
CB07 model with the bluest colors in the figure. In general, 
the objects in the z ~ 8 sample have fainter magnitudes 
than the objects in the z ~ 7 sample. Along with the smaller 
sample size, the measurements of the mean colors of the z ~ 

8 sample have larger uncertainties. 

4.2. UV Spectral Slope 

In order to facilitate the comparison of our results to pre- 
vious studies, we use the WFC3 photometry to approximate 
the slope of the rest-frame UV continuum, /?, of the galax- 
ies in our z ~ 7 sample, where f\ cx A' 3 (e.g., Calzetti et al. 
1994). Although the definition of (5 is based on regions of 
the UV spectrum devoid of strong absorption features, it is 
common to estimate the UV spectral slope of distant galaxies 
using the broad-band measures of the rest-frame UV color 
(e.g. Meureret al. 1999). Here we use the relation of Bouwens 
et al. (2009b) to estimate (3 for galaxies in our z ~ 7 sample, 

/3p ho t = 4.29(/i 2 5-#i6o)-2.0, (3) 

where /3 p h t denotes a photometric estimate of /?. The formal 
definition of j3 pertains only to wavelengths greater than that 
of Lya. Therefore, we restrict our analysis to those objects 
in our z ~ 7 sample with z p h t < 7.5, where Lya emission 
falls outside the F125W-band. We show /3 p h t for our sample 
at Zphot < 7.5 in Figure 6 along the middle vertical axis, and 
our measurements of /3 p hot are tabulated in Table 2, from the 
colors from the Monte Carlo bootstrap simulations. 

4.3. Interpretation of the UV Colors and UV Slopes 

Figure 6 compares our measurements of the mean color and 
its associated uncertainty for the z ~ 7 sample to the empiri- 
cal Kinney et al. (1996) local starburst template spectra. We 
also compare our measurements to stellar-population synthe- 
sis models from CB07 spanning a range of age, metallicity 
and dust extinction. Lastly, we also compare our measure- 
ments to the local extremely blue galaxy NGC 1705, also 
from Kinney et al. (1996). NGC 1705 is a blue compact dwarf 
galaxy with recent strong star-formation activity, and is unob- 
scured with E(B-V) = 0.00 (Calzetti et al. 1994). Lee & Skill- 
man (2004) published metallicity measurements from 16 H II 
regions in NGC 1705, finding an average oxygen abundance 
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FIG. 6. — //i60 versus J125 -ffieo color-magnitude diagram for the z ~ 7 sample. The 6.3 < z p h t < 7.5 galaxies are plotted as gray squares, with the error 
bars representing the 1 a photometric uncertainties. The gray shaded area corresponds to the region excluded by the magnitude limits of our sample. The central 
vertical axis denotes the UV spectral slope f3, which is a function of the color (see §4.2). The black, red and blue large circles show the mean color for all 
objects, the bright subsample (//160 < 28.3) and the faint subsample (//|60 > 28.3), respectively. The dotted line denotes the magnitude split of H\(,o = 28.3. 
The histograms in the right-hand panel represent the distribution of the means for all objects (black histogram), objects with H^q < 28.3 (red histogram), and 
objects with //jgo > 28.3 (blue histogram) from the Monte Carlo bootstrap simulations. The errors on the means were derived from the standard deviation of 
these distributions. The gray cross-hatched region shows the color spanned by the Kinney et al. (1996) starburst galaxy templates (redshifted to z = 7). We also 
show the Kinney template for NGC 1705, which is one of the bluest local starburst galaxies. The colored arrows denote the colors of CB07 stellar population 
models, described in the figure legend, where the three numbers denote their age, metallicity (in solar units), and dust extinction (in Ay), respectively. While 
many objects have measured colors bluer than NGC 1705, the average of the faint subsample is consistent within 2 a. The faint subsample is also consistent with 
all dust-free CB07 models. This implies that these objects may have very low dust contents. However, as the distribution of the CB07 models shows, the colors 
of stellar populations do not change much between 1-100 Myr and Z = 0.005-1.0 Zq, so determinations of these properties are difficult, though it does appear 
as if these objects are < 100 Myr and/or have Z < Zq. 



of ~ 0.35 Zq, making it comparable to Local Group dwarf 
irregular galaxies. More recently, Annibali et al. (2009) find 
a good fit to the young stellar populations in NGC 1705 when 
assuming metallicities of only 0.02 Zq. Thus, NGC 1705 is 
one of the more metal poor local star-forming galaxies. 

Both the bright and faint z ~ 7 subsamples are bluer than 
local starburst templates from Kinney et al. (1996). With our 
simulation results, we derive a 99.9987% (4.4 er) confidence 
that the color of an average z ~ 7 galaxy is bluer than that of 
a local starburst galaxy. The confidence that the bright and 
faint subsamples are individually bluer than local starbursts 
is slightly less given their smaller sample sizes, at 99.88 % 
and 99.86 % (~ 3.2 a), respectively. Rather, the colors of 
these z ~ 7 galaxies are closer to the that of NGC 1705. The 
confidences that these three samples are bluer than NGC 1705 
are ~ 80%, 33% and 95%, respectively (i.e. consistent within 

2(7). 

This is in contrast to studies of star-forming LBGs at lower 
redshifts, 2 < z < 6 (e.g., Papovich et al. 2001; Ouchi et al. 
2004; Reddy & Steidel 2009; Overzier et al. 2009; Bouwens 
et al. 2009a). Papovich et al. (2001) showed that UV-selected 
galaxies at z ~ 2-3 have rest-frame UV colors which span 
a range similar to that of the composite starbursts from Kin- 
ney et al., while these typical local templates are unable to 
reproduce the blue colors of the z ~ 7 sample. Therefore, we 
conclude that z ~ 7 galaxies are dominated by stellar popula- 
tions that have a combination of young ages, lower metallic- 



ity, and little dust extinction. Additionally, in contrast to the 
color-magnitude relations observed in lower-redshift galax- 
ies at 4 < z < 6 (e.g., Papovich et al. 2004; Bouwens et al. 
2009a; Overzier et al. 2009), with our relatively small sample 
size we observe only a weak (< 2 a) dependence between the 
UV-color and magnitude in the z ~ 7 sample. 

Galaxies at z ~ 7 are uniformly blue, nearly independent of 
magnitude. Figure 6 shows the expected colors of six CB07 
stellar population synthesis models, covering a range in age 
(1-100 Myr), metallicity (0.005-1 Z Q ), and dust extinction 
(Ay = 0-1 mag) formed with a constant star-formation his- 
tory. With only a single rest-frame UV color, it can be dif- 
ficult to accurately constrain these properties. For example, 
for a fixed age of 1 Myr (100 Myr), changing the metallicity 
from 0.005 Zq to Z© increases the ./125-H160 by 0.08 mag 
(0.16 mag). Similarly, at fixed metallicity increasing the age 
from 1 Myr to 100 Myr increases the color by 0.07-0. 14 mag, 
with a slight dependence on the metallicity. The models with 
either very young ages (< 100 Myr) or very low metallic- 
ity (~ 0.005 Zq) match the observed /125 -#160 colors of the 
faint z ~ 7 subsample with > 28.3 mag. This is not true 
for models with higher ages or higher metallicities. If we in- 
clude the effects of dust attenuation on any of these models 
the colors become considerably redder, extending them be- 
yond the la range of the faint Hi^o > 28.3 mag subsample, 
and also extending them beyond the 1 <7 range for the brighter 
Hi 60 < 28.3 mag for the model with age 100 Myr. The conclu- 
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FIG. 7. — H[(,o versus J\2s 
and models redshifted to z = 



//l60 color-magnitude diagram for the z ~ 8 sample. The Zph ot > 7.5 galaxy candidates are plotted with symbols defined in Figure 6, 
!. As hya affects the F125W-band at z ~ 8, we show an additional model track as a dashed line, which is the bluest possible model 
from CB07, with hya added in (this model has a rest-frame hya EW ~ 430 A; see §5.1). Similar to z ~ 7, these objects appear to be fairly dust free, especially 
the faint subsample, but determinations of their metallicities and ages are difficult, as the object's colors are consistent with a range of both parameters. At these 
redshifts hya lies in the F125W-band, thus we do not compute /3. 



sion here is strong. The extremely blue colors of faint galaxies 
with Hieo > 28.3 mag in the z ~ 7 sample implies that they 
have low dust contents. Their colors are also consistent with 
low metallicities and low stellar population ages, but further 
work is necessitated to place strong constraints on these prop- 
erties (see §5). 

We draw similar conclusions for the z ~ 8 sample as for 
the z ~ 7 sample. However, due to the smaller sample size 
and increased distance to these galaxies, the conclusions have 
a lower significance, as the bin of all galaxies are bluer than 
local starbursts at only ~ 1.7 sigma significance. The mean 
color J\ 25- H\ 60 =-0.01 ±0.22 mag of the faint subsample is 
consistent with that at z ~ 7 within ~ 1 er, and it is also con- 
sistent with the expected color of NGC 1705. 

Schaerer (2002) have argued that nebular continuum emis- 
sion from regions around massive star formation at very 
young ages could affect the observed colors of high-redshift 
galaxies at these redshifts. While this continuum emission is 
not included in the CB07 stellar population synthesis models, 
we conclude here that its effects are mitigated in our sample 
(similar to conclusions reached by Bouwens et al. (2009b)). 
The local starburst galaxy I Zw 18 has the lowest metallicity 
known of any galaxy, and it is dominated by actively forming 
young stars as evidenced by its large Wolf-Rayet stellar pop- 
ulation. Nevertheless, I Zw 18 shows no evidence for nebular 
continuum emission (Brown et al. 2002). Furthermore, while 
this continuum emission is expected to redden the UV con- 
tinuum of the galaxies, mimicking stellar populations with a 
combination of higher dust attenuation, older ages, or higher 
metallicity, the galaxies in our z ~ 7 and z ~ 8 samples have 
nearly maximally blue colors, implying that this red contin- 
uum does not dominate their integrated properties. Lastly, 
the nebular continuum is significant only if a high fraction of 
hydrogen-ionizing photons are absorbed by the ISM, which 



may not be the case here as the dust extinction appears small 
(see also Bouwens et al. 2009b, and § 7). This may imply 
that the fainter galaxies have higher Lyman continuum escape 
fractions (see § 7). 

4.4. Comparison with Other Samples 

Several groups have studied the slope of the rest-frame 
UV continuum of LBGs (e.g., Steidel et al. 1999; Meurer 
et al. 1999; Ouchi et al. 2004; Stanway et al. 2005; Reddy 
et al. 2006, 2008; Reddy & Steidel 2009; Hathi et al. 2008; 
Bouwens et al. 2006, 2007, 2009a), finding a consensus con- 
clusion that the UV spectral slope evolves toward lower val- 
ues (bluer colors) at progressively higher redshifts and fainter 
magnitudes. Bouwens et al. (2009a) found that LBGs at at 2.5 
< z < 6 shift systematically to bluer colors with increasing 
redshift, even when properly taking into account the LBG- 
selection criteria. The slope of the U V continuum (3 should be 
more sensitive to variations in the amount and configuration 
of dust attenuation in these galaxies, and rather less sensitive 
to differences in stellar ages, metal fraction, or the IMF (e.g., 
Calzetti et al. 1994; Meurer et al. 1999). Therefore, this evolu- 
tion in (3 implies that the overall dust content in distant galax- 
ies decreases with increasing redshift. Additionally, Bouwens 
et al. (2009b) reported evidence for a shift of /3 p h t to steeper 
(i.e., bluer) values at lower UV luminosity, as high-redshift 
galaxies with L ~ 0.15 L* 3 are systematically bluer by ~ 0.3 
in Pphot compared to L! =3 galaxies, where L! =3 refers to the 
characteristic luminosity at z = 3 measured by Steidel et al. 
(1999). 

For our full sample of galaxies with z < 7.5 we find a mean 
UV spectral slope of /3 p h ot = -2.44 ± 0.25, dropping to /3 p h ot = 
-3.07 ± 0.51 for the faint subsample with Hk,q > 28.3 mag. 
Evidence for a dependence of UV slope on the UV luminosity 
is only weakly observed in our data given the uncertainties. 
These values for the spectral slopes are indicative of galax- 
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TABLE 2 

Average UV Continuum Properties 



Redshift 


Sample 


7l25 -#160 


/3phot 


Confidence bluer than 


Confidence bluer than 


Puv 








local starbursts 


NGC1705 


(x 10 26 ) 



Z r 


- 7 


All Objects 


-0.10 ±0.06 


-2.44 ± 0.25 


99.9987% 


79.8% 


0.44 ±0.11 


Z r 


7 


#160 < 28.5 


-0.03 ± 0.06 


-2.13 ±0.25 


99.88% 


32.7% 




Z r 


-7 


#160 > 28.5 


-0.25 ±0.12 


-3.07 ±0.51 


99.86% 


95.0% 




Z r 


- 8 


All Objects 


0.00 ±0.16 




91.1% 


20.3% 


0.14 ±0.06 


Z r 


- 8 


#160 < 28.5 


0.08 ±0.19 




74.9% 


14.2% 




Z r 


- 8 


#160 > 28.5 


-0.01 ± 0.22 




84.6% 


29.9% 





NOTE. — We do not measure /3 p hot for the z ~ 8 galaxies, as Lya (emission or absorption) will contaminate 
the measurement. Additionally, the confidence ranges on the z ~ 7 galaxy sample are more robust than at z ~ 8, 
given the larger number of objects in the z ~ 7 sample. We derive the errors on these values via 10 7 bootstrap 
Monte Carlo simulations. The specific luminosity density at 1500 A, puv, has units of erg s -1 Hz -1 Mpc~ 3 . 



ies dominated by young stars with low dust attenuation (e.g., 
Calzetti et al. 1994; Kinney et al. 1996), and thus reinforces 
our conclusions from the §4.3. In particular, the low values 
of /3 p hot for our faint subsample are consistent with the local 
starburst NGC 1705, which has near-zero dust attenuation. 

Bouwens et al. (2009b) obtain similar UV spectral slopes in 
their study of the UV continuum slopes of their z ~ 7 sample 
(with -20 < M1500 < -18). Bouwens et al. argue that this is 
evidence for very low metallicities (Z < 0.005 Zq, possibly 
as low as Z = 0) as only these stellar populations reproduce 
the observed values for j3. The models studied by Bouwens et 
al. include effects of nebular continuum (see discussion in the 
previous section), which would redden the UV slope, imply- 
ing that the galaxies would need a very large escape fraction 
for UV photons (/ ejc > 0.3). However, as we discuss in §4.3, 
the distribution of galaxy colors after properly accounting for 
the photometric uncertainties does not provide evidence for 
this effect. Rather, the evidence suggests these objects have 
a make-up consistent with extremely blue local galaxies such 
as NGC 1705. 

5. THE STELLAR POPULATIONS OF 6.3 <z< 8.6 GALAXIES 
5.1. Stellar-Population Synthesis Model Fitting 

We compare the galaxy colors to stellar population synthe- 
sis models to further constrain the properties of these galaxies. 
We fit all 3 1 galaxies in our z ~ 7 and z ~ 8 samples with a 
suite of CB07 models, and in all cases we use models with 
a Salpeter IMF (with the default CB07 mass range of 0.1 - 
100 Mq). Because the rest-frame UV colors of the galaxies 
are indicative of ongoing star formation we have fixed the star 
formation history (SFH) to be constant (see also Labbe et al. 
2009a). Given the young ages we derive for the stellar pop- 
ulations, this is a reasonable approximation. However, if we 
allow for exponentially declining SFHs, we would expect to 
find even younger stellar population ages (see, e.g., Papovich 
et al. 2001; Shapley et al. 2001), therefore this assumption 
yields conservative results. 

During the fitting, we vary the age (where age represents 
the time elapsed since the onset of star formation) from 1 Myr 
to the time elapsed from z = 20 (a reasonable estimate for the 
onset of galaxy formation) to the given redshift, and we apply 
dust extinction using the Calzetti et al. (2000) attenuation law 
over a range of < E(B- V) < 0.4 (equivalent to < A v < 1 .6 
mag). Given the blue colors of these objects from §4, we ex- 
pect our maximum possible dust extinction to be sufficiently 



high. We allow the metallicity to vary between the full range 
in CB07, from 0.005 - 2.5 Z . We apply IGM attenuation via 
the prescription of Madau (1995). 

We expect the galaxies in our samples to be dominated by 
intensely star-forming stellar populations. We therefore in- 
clude the effects of Lya emission on the stellar population 
models (for details see Finkelstein et al. 2008, 2009). Lya 
emission may significantly affect the colors given that the ob- 
served EW of the line increases as (1 + z) and the effective 
width of the filters decreases as (1 +z) _1 . To estimate the ex- 
pected Lya line-flux, we use the number of ionizing photons 
provided for each CB07 model as a function of stellar popula- 
tion age. Assuming Case B recombination (Osterbock 1989), 
approximately 2/3rds of hydrogen-ionizing photons produce 
Lya photons during recombination. We sum the expected 
Lya flux with the CB07 model spectrum. This is done prior 
to the application of the dust and IGM attenuation in order 
to most accurately represent the processes in these galaxies. 
Lya lines in Lya emitting galaxies typically show line profiles 
truncated on the blue-side presumably due to IGM absorption 
(e.g., Stern & Spinrad 1999; Rhoads et al. 2003). We model 
this effect by allowing the full IGM attenuation on one-half 
of the line flux. This results in an effective IGM optical depth 
at Lya of 

Teff = -ln(0.5)-ln(l +exp[TiGM(Az, V a)]). (4) 

We tested the redshift dependence on our SED fits by al- 
lowing the redshift to be a free parameter, fitting the full ACS, 
WFC3 and IRAC photometry over a grid of < z < 11, with 
6z = 0.05. We find that differences in stellar libraries and IGM 
treatment between our fitting method and that of EAZY re- 
sult in a slightly higher redshift from this fitting, by a few 
tenths on average, with a range in the sample from 6.50 < z 
< 8.85. At these correspondingly larger distances, the best- 
fit models are brighter, which the fitting process compensates 
for by making the models younger (in many cases maximally 
young; 1 Myr). Younger stars are more luminous in the UV, 
thus the corresponding stellar mass is less. When compar- 
ing the 68% confidence ranges from this fitted redshift to the 
photometric redshift, we find that they are consistent in nearly 
all cases. We thus adopt z = z p hot for the model fitting, as this 
will again yield conservative results (i.e., older and more mas- 
sive). However, we will examine how the uncertainties on this 
redshift affect our parameter uncertainties below. 

We fit the model spectra, integrated with the photometric 
bandpasses (including bandpass transmission, telescope op- 
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tics and detector efficiency) to the photometry from zgso, T105, 
J\25, #160, and IRAC [3.6] and [4.5] for each galaxy. For each 
model we calculate x 2 using the photometric errors, including 
a systematic error of 0.05 x/„, where f v is the flux measured 
in each bandpass. This systematic error accounts for possible 
zeropoint uncertainties, aperture corrections, and other uncer- 
tainties that scale multiplicatively with the flux. 

During the fitting, any data point with a wavelength below 
Lya detected with < 2 a significance was not included in the 
fit, unless a given model violated the 2 a upper limit, in which 
case the \ 2 for that band is added in, using the 2 a upper limit 
as the flux, and the 1 a flux as the error. At the photometric 
redshift, the flux in all models should be heavily attenuated 
by the IGM, so this essentially results in these points being 
excluded by the fit, which is preferable as there can be large 
variations in the IGM optical depth along various lines-of- 
sight. 

For the IRAC data most galaxies are undetected, and we al- 
lowed the models to fit the actual measured fluxes, as the pho- 
tometric errors still provide constraints on the possible model 
parameters. However, their contribution to the % 2 is minimal 
due to their large uncertainties, thus the IRAC fluxes for these 
galaxies for the most part do not provide strong constraints on 
the shape of the SED. Rather, the blue rest-frame UV (Y105 - 
7i25 and J\25 - #160) colors coupled with the constraint that the 
stellar population age be less than the time elapsed since z = 
20 places stronger constraints on the stellar population prop- 
erties. 

Although these objects are very blue in the rest-frame 
UV, we are unable to exclude the possibility that the galax- 
ies contain generations of stars from previous star-formation 
episodes, lost in the "glare" of the UV-luminous popula- 
tion. We allowed for the possibility that we were missing an 
older generation of stars by generating a second set of stellar- 
population synthesis models with two components. The first 
model component is similar as described above, and is used 
to describe the UV emission from young stars, except that 
the star formation history is modeled as an instantaneous star- 
burst. The second stellar population is also modeled as an 
instantaneous starburst, with an age equal to the time elapsed 
between z = 20 and z p h t derived for each galaxy. Because of 
the limited number of data points we have for each galaxy, we 
considered only models where a fixed mass fraction of 90% 
of the mass exists in this "maximally old" stellar population. 

5.2. Fitting Results 

Figures 8 and 9 display the best-fit models for the galax- 
ies in our z ~ 7 and z ~ 8 samples, respectively, and Table 3 
provides tabulated results. The best-fit model parameters re- 
fer to those models with the minimum \ 2 calculated between 
the object and model fluxes. We calculate the 68% confidence 
ranges on each of the model parameters (age, dust, metallicity, 
and stellar mass) for our fits using a series of 10 3 Monte Carlo 
bootstrap simulations including the photometric uncertainties 
and the probability distribution function of the photometric 
redshifts. The resulting confidence range is the central 68% 
of the best-fit results for each parameter for each object. In 
each simulation we vary the input fluxes by an amount pro- 
portional to their errors, such that the parameter uncertainties 
reflect the photometric uncertainties. Table 3 lists the 68% 
confidence range on the parameters for each object. The best- 
fit ages of these objects from the full sample range from 1 - 
720 Myr, with a median of ~ 200 Myr. The best-fit color ex- 



cesses range from < A v < 1.5 mag (0 < E(B-V) < 0.36), 
with a median of Ay — 0.3. For individual galaxies there is 
a relatively large 68% confidence range on the model stellar- 
population age, dust color excess, and metallicity. This is due 
to the fact that these parameters are degenerate (as they all 
result in reddening of the spectra), and the large 68% range 
reflects this degeneracy. In the absence of data measuring the 
strength of Balmer/4000 A-break, we are unable to improve 
these constraints. In contrast, our model fitting provides rela- 
tively tighter constraints on the galaxy stellar masses. 

To quantify the effect that the redshift uncertainties have on 
our derived parameter uncertainties, we also performed a set 
of simulations assuming a redshift uncertainty of zero. From 
both sets of simulations, we compute characteristic parame- 
ter errors by taking the mean of the difference between the 
best-fit value and the lower and upper 68% confidence range. 
When we compare the stellar population parameter uncertain- 
ties with and without redshift variation, we find that account- 
ing for the redshift uncertainty increases the errors by ~ 30% 
in mass, 10% in age, and 25% in E(B-V). We thus conclude 
that the effect of redshift uncertainties should be included in 
the error propogation when using photometrically estimated 
redshifts. In the subsections that follow we discuss the con- 
straints on the various stellar population parameters. 

5.2.1. Stellar Population Age, Dust Extinction and Metallicity 

The median best-fit age in our sample is ~ 200 Myr, but all 
objects have 68% confidence ranges which are consistent with 
ages < 80 Myr, with 24/3 1 objects consistent with extremely 
young ages of < 10 Myr. However, the large spread in ages 
make it difficult to make strong conclusions. Nevertheless, the 
colors of the majority of galaxies are consistent with model 
fits requiring young stellar-population ages. 

The best-fit dust extinctions for galaxies in our sample 
range from < A v < 1.5 mag, with a median extinction of 
A v ~ 0.3 mag. However, the 68% confidence range on the 
extinction for 16/31 objects are consistent with zero dust ex- 
tinction. This is consistent with our conclusions in § 4 based 
on the very blue colors of these objects. 

The fact that a few galaxies show redder UV colors implies 
some dust attenuation, especially at brighter //160 magnitudes. 
This is consistent with dust-production rates in core-collapse 
supernovae (SNe). Todini & Ferrara (2001) show that ~ 0.1- 
0.3 M of dust form per star with initial masses of 12-35 
M Q , which end as SNe after < 20 Myr. Taking a typical 
stellar mass of 10 9 M Q for galaxies in our z ~ 7 and z ~ 8 
samples (see § 5), we expect roughly ~ 2 x 10 6 stars in this 
mass range will form (assuming a Salpeter IMF). This implies 
the production of nearly 10 6 M© of dust from core-collapse 
supernovae. Assuming this dust mixes with the galaxies ISM 
on a dynamical timescale (~ 30 Myr assuming a size of 1 .5 
kpc and a circular velocity of 50 km s" 1 ), it is highly plausible 
that this dust has a measurable effect on the z ~ 7 and z ~ 8 
galaxy rest-frame UV colors. 

Surprisingly, our model fits place strong constraints on the 
metallicities of the galaxies in our samples, relative to the 
age and extinction. The majority (19/31) of galaxies in our 
samples have best-fit metallicities Z = 0.005 Z©. A further 
8 of the remaining 12 galaxies have best-fit metallicities of 
Z = 0.02 Z©. Additionally, all 31 galaxies have their entire 
68% confidence ranges at Z < 0.02 Z Q . This is consistent 
with our conclusions in § 4 that the metallicites of the major- 
ity of galaxies at z ~ 7 and 8 are a fraction of solar. Galaxies 
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FIG. 8. — The observed spectral energy distributions (SEDs) of our objects at 6.3 < Zphot < 7.5 with the best-fit model spectra overplotted, shown in order of 
increasing Zphot- The spike in the spectra represents Lya emission as described in the text. The best-fit results are tabulated in Table 3. When a data point is 
detected at less than 2 a significance, we show the 1 cr upper limits as arrows. The blue colors of these objects typically result in best-fit spectra that are blue, 
implying very young ages. We show the best-fit stellar mass in the top of each panel. 
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FIG. 8. — Continued. 
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TABLE 3 
Model Fitting Results 



ID Zp), ot Mass Mass Max Mass Max Mass Age Age Ay Ay Z Z 

Best-Fit 68% Range Best-Fit 68% Range Best-Fit 68% Range Best-Fit 68% Range Best-Fit 68% Range 
(1O 8 M ) (10 s M ) (10 s M ) (10 s M ) (Myr) (Myr) (mag) (mag) (Z ) (Z ) 



803 


6.34 


14.2 


3.9- 


25.6 


20.7 


4.6 


-26.1 


700 


50 


-700 


0.32 


0.08 


-0.89 


0.005 


0.005 


- 0.200 


2894 


6.34 


8.9 


2.6- 


15.2 


15.2 


3.9 


- 16.6 


400 


50 


-700 


0.24 


0.00 


-0.65 


0.005 


0.005 


- 0.020 


769 


6.36 


18.4 


12.1 - 


44.8 


25.0 


18.7 


-74.5 


100 


10 


- 100 


0.49 


0.49 


- 1.38 


2.500 


0.020 


-2.500 


1289 


6.38 


6.1 


1.2- 


12.2 


5.9 


2.3 


- 12.4 


300 


30 


-700 


0.24 


0.00 


-0.65 


0.005 


0.005 


- 0.400 


649 


6.40 


5.6 


2.9- 


11.9 


24.3 


9.8 


-33.4 


60 


20 


-200 


0.65 


0.32 


-0.89 


0.005 


0.005 


- 0.020 


1445 


6.40 


3.2 


0.3- 


5.7 


5.0 


2.2 


- 6.5 


400 


5 


-500 


0.00 


0.00 


-0.40 


0.005 


0.005 


- 0.005 


2032 


6.40 


7.9 


0.7- 


12.3 


5.4 


2.3 


- 10.7 


600 


5 


-700 


0.08 


0.00 


-0.57 


0.005 


0.005 


- 0.200 


1818 


6.43 


0.2 


0.2- 


3.6 


2.5 


2.0 


- 8.0 


1 


1 


- 60 


0.00 


0.00 


-0.81 


2.500 


0.005 


- 2.500 


328 


6.47 


3.7 


0.2- 


6.0 


2.9 


1.2 


-11.8 


700 


1 


-700 


0.00 


0.00 


-0.97 


0.020 


0.005 


- 0.200 


1072 


6.47 


7.1 


0.3- 


7.7 


5.1 


2.1 


- 8.9 


700 


2 


-700 


0.08 


0.00 


-0.65 


0.005 


0.005 


- 0.200 


1566 


6.56 


20.8 


10.1 - 


31.8 


37.8 


16.8 


-41.6 


300 


80 


-600 


0.49 


0.24 


- 0.73 


0.005 


0.005 


- 0.005 


2013 


6.70 


2.6 


0.4- 


6.0 


5.9 


2.8 


- 10.1 


200 


8 


-400 


0.08 


0.00 


-0.57 


0.005 


0.005 


- 0.005 


1110 


6.72 


4.5 


0.6- 


10.8 


8.4 


3.2 


-31.1 


20 


1 


- 30 


0.97 


0.16 


- 1.38 


0.020 


0.005 


-2.500 


1441 


6.83 


17.5 


10.5- 


29.4 


32.7 


20.8 


-43.2 


20 


20 


- 40 


0.73 


0.57 


-0.89 


0.020 


0.020 


- 0.020 


567 


6.84 


0.3 


0.2- 


3.0 


1.8 


1.7 


- 8.4 


10 


1 


- 30 


0.16 


0.00 


-0.89 


0.020 


0.005 


- 2.500 


2432 


6.86 


7.7 


2.4- 


12.6 


22.9 


5.8 


-26.9 


60 


8 


-100 


0.57 


0.24 


-0.81 


0.005 


0.005 


- 0.400 


515 


6.88 


0.6 


0.2- 


3.7 


1.4 


1.7 


- 6.2 


60 


1 


-100 


0.00 


0.00 


-0.73 


0.005 


0.005 


- 0.200 


335 


6.93 


6.3 


0.6- 


10.6 


10.0 


2.4 


- 18.6 


400 


4 


-600 


0.32 


0.00 


- 1.05 


0.005 


0.005 


- 0.020 


1768 


7.22 


8.1 


2.4- 


14.8 


9.9 


9.5 


-33.4 


10 


3 


- 30 


0.97 


0.57 


- 1.05 


0.020 


0.020 


- 2.500 


2056 


7.25 


20.6 


17.3- 


95.6 


85.5 


53.8- 


-135.3 


5 


4 


-600 


1.46 


0.65 


-1.62 


0.005 


0.005 


- 0.200 


669 


7.25 


5.4 


2.5- 


24.9 


11.7 


7.8 


-52.8 


7 


1 


-200 


1.05 


0.49 


- 1.46 


0.020 


0.005 


-2.500 


1106 


7.32 


7.8 


1.8- 


14.7 


14.7 


6.8 


-22.0 


200 


20 


-400 


0.16 


0.00 


-0.57 


0.005 


0.005 


- 0.020 


3053 


7.40 


5.0 


0.7- 


11.1 


9.4 


3.8 


- 16.7 


200 


S 


-500 


0.16 


0.00 


-0.81 


0.005 


0.005 


- 0.020 


819 


7.55 


16.2 


1.5- 


21.8 


13.1 


4.3 


-36.7 


500 


2 


-500 


0.40 


0.16 


- 1.30 


0.005 


0.005 


- 0.200 


653 


7.76 


0.2 


0.2- 


2.4 


2.0 


1.9 


- 2.7 


4 


1 


-200 


0.00 


0.00 


-0.24 


0.005 


0.005 


- 0.200 


2055 


7.81 


8.8 


6.6- 


31.1 


43.8 


25.2 


-88.8 


2 


2 


- 10 


1.38 


1.21 


- 1.62 


2.500 


0.020 


- 2.500 


200 


7.91 


6.4 


0.3- 


9.4 


3.5 


2.8 


- 13.7 


500 


1 


-500 


0.00 


0.00 


-0.81 


0.005 


0.005 


- 0.200 


213 


8.05 


5.3 


0.2- 


14.7 


10.5 


1.8 


-31.5 


400 


1 


-400 


0.16 


0.00 


- 1.54 


0.200 


0.005 


- 2.500 


3022 


8.05 


19.8 


0.3- 


50.9 


20.6 


1.7 


-66.1 


400 


1 


-400 


0.89 


0.16 


-1.62 


0.020 


0.005 


- 2.500 


640 


8.24 


6.8 


1.4- 


12.8 


18.8 


7.4 


-46.6 


20 


1 


- 20 


0.65 


0.16 


- 1.13 


0.020 


0.005 


- 2.500 


125 


8.61 


4.3 


0.8- 


20.6 


9.1 


3.5 


-43.0 


40 


2 


-300 


0.73 


0.24 


- 1.62 


0.005 


0.005 


-2.500 



NOTE. — We derive the 68% confidence range on our parameters by taking the central 68% of the probability distribution for 
each parameter from the 10 3 bootstrap Monte Carlo simulations, including the photometric uncertainties and the photometric redshift 
probability distribution function. The maximum mass is the best-fit mass from the two-component models. 



at these redshifts do not have UV colors consistent with so- 
lar metallicities. Figure 10 shows the distribution of best-fit 
metallicities for our samples, as well as the joint probabil- 
ity distribution function for metallicity, constructed from the 
probability distribution function for each object. Only galax- 
ies with brighter apparent magnitudes, H\&q < 28.3 mag have 
possibly solar (or super-solar) best-fit metallicities. Given 
the model-parameter uncertainties, the joint probability dis- 
tribution is more informative than the individual best-fits. The 
probability distribution has a maximum at Z = 0.005 Z© with a 
wide distribution spanning the full parameter space. Figure 10 
also shows the cumulative probability distribution function. 
We find that 68% of the probability distribution function is 
constrained at Z < 0.05 Z Q . Furthermore, 95% of the cumu- 
lative probability distribution function lies within metallicities 
< Zq . Therefore, we conclude that the majority of galaxies in 
our sample are best-fit by sub-solar metallicities, and in some 
cases are only a few percent of solar. 

5.2.2. Stellar Mass 

The stellar masses for each object are typically better con- 
strained compared to the other model parameters. In general 
these constraints result from the fact that the relatively young 
universe (< 800 Myr) limits the age of the stellar populations, 
and thus limits the amount of stellar mass in "evolved" stellar 
populations. With this constraint, the rest-frame UV colors 
alone are sufficient to provide constraints on the total mass- 



to-light ratios and therefore the stellar masses. 

As given in Table 3, the median 68% confidence range on 
stellar mass spans factors of ~ +2/— 4.5. Thus, the upper limit 
on the masses is constrained to within a factor of two, while 
they may be a factor of 4.5 less massive. Figure 1 1 shows 
the distribution of best-fit stellar masses for the z ~ 7 and z ~ 
8 populations. We find average masses of slightly less than 
~ 10 9 Mq for galaxies at z ~ 7 and z ~ 8, consistent with the 
stacking analysis of z ~ 7 galaxies from Labbe et al. (2009b). 
The figures also show the joint probability distribution for the 
stellar masses, constructed using the results for each individ- 
ual galaxy. The peak in the joint probability distributions is 
roughly the same at z ~ 7 and z ~ 8, with a maximum at 
« 7 x 10 8 M Q , though z ~ 8 LBGs are more likely to have M 

- 10 7 Mq. 

5.2.3. Maximal Stellar Mass 

Figure 1 1 also shows the distributions for the "maximal" 
stellar mass, defined as the stellar mass of our two-component 
models where 90% of the stellar mass exists in a passively 
evolving population formed instantaneously at z = 20. Sur- 
prisingly, we find that the maximal mass is shifted to a median 
value of slightly less than twice that for the best-fit single- 
population model. Additionally, in many cases the upper end 
of the 68% confidence range on the best-fit single-population 
model masses exceed the maximum mass value. However, the 
most massive of either quantity is < 10 10 M Q and they have 
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FIG. 10. — The histogram shows the best-fit metallicites from our model 
fitting, with red and blue lined regions denoting //i60 < 28.3 and //i60 > 
28.3, respectively. The black line denotes the metallicity joint probability 
distribution for the whole sample. The dashed line shows the cumulative 
probability, integrating from the peak in the probability distribution at Z = 
0.005 Zq . This curve changes color from blue to green to red while in the 
68, 95 and 99.5% confidence ranges, respectively. From this, we rule out 
solar metallicities at 95% confidence for the sample, and constrain Z < 0.05 
Zq with 68% confidence (c.f., Bouwens et al. 2009b). 

similar distributions, thus in the following discussion, we use 
the best-fit two-burst mass as an indicator of the maximum 
mass. This distribution is shifted only marginally to higher 
masses compared to the mass distribution from the single- 
component models discussed in the previous section. This 
small shift is counterintuitive as our constraints stem primar- 
ily from photometry of the rest-frame UV, and contrasts with 
studies of lower redshift LBGs, which have maximal masses 
larger by up to an order of magnitude (e.g., Papovich et al. 
2001; Erb et al. 2006b; Pirzkal et al. 2007; Finkelstein et al. 
2009). 

To increase the maximal mass would require an older stellar 
population with a higher mass-to-light ratio. However, this 
would require ages greater than the lookback time from z = 
20 to the observed redshift. Allowing the stellar mass in the 
older stellar population to increase would increase the emitted 
flux from this population in the F125W and F160W bands, 
increasing the \ 2 an d reducing the quality of fit. To verify 
this, we also performed bootstrap Monte Carlo simulations 
on the two-burst models. If the flux contributed from the old 
burst to the WFC3 bands was negligible, one would expect 
the 68% confidence range on the maximal mass to be quite 
large. However, as shown in Table 3, we find this to not be the 
case, as the maximal masses are quite well constrained, with 
a median uncertainty of a factor of ± 2.5 

If the stellar mass fraction in old stars was above 90%, this 
would also increase the total stellar mass. Allowing 99% of 
the stellar mass to be in this component, the maximal stellar 
masses increase by a median factor of ~ 3, not a factor of 10 
as may be expected if the choice of mass ratio was driving our 
results. However, we consider these masses unlikely as our 
model already adopts an extreme case where the vast majority 
of stellar mass formed at very high redshift. If this is true for 
every galaxy, then we would expect a significant population 
of galaxies at z ~ 7 and z ~ 8 dominated by > 500 Myr stellar 
populations, which is not present. Given the young ages we 
derive for the z ~ 7 and z ~ 8 samples (a few x 10 Myr) 
there would need to be a very high duty cycle of rejuvenated 
star formation every few x 10 Myr to reproduce the current 
observations. 



We thus conclude that at z ~ 7 and z ~ 8, the relatively 
young age of the Universe since z = 20 (400-700 Myr) con- 
strains the age of any pre-existing stellar population. This lim- 
its the contribution to the mass-to-light ratio to come solely 
from from earlier-type stars (earlier than A-type stars with 
> 3 M ), and limits the amplitude of the 4000 A break. 
The IRAC data available for the HUDF provides generally 
a weaker constraint on the stellar mass in these galaxies rela- 
tive to the constraint from the age of the Universe (c.f., Labbe 
et al. 2009b). 

5.3. Comparison of Stellar Masses to the Literature 

The interquartile range of stellar masses (encompassing the 
inner 50%-tile) for our z ~ 7 and 8 samples spans 4.5-8.9 x 
10 8 M Q . These are lower than the stellar masses inferred for 
UV-selected galaxies at 2 < z < 6. At z ~ 2- 3, the inferred 
stellar mass of "L*" UV-selected galaxies are on the order of 
10 10 M Q (Papovich et al. 2001; Shapley et al. 2001; Erb et al. 
2006b; Reddy et al. 2006), and there are UV-selected galax- 
ies with stellar masses exceeding 10 11 Mq at these redshifts 
(Shapley et al. 2005; Erb et al. 2006b; Forster Schreiber et al. 
2009). The fact that the galaxies at 6.3 < z < 8.6 in our sample 
have such low masses implies a rapid decrease in the stellar 
mass density of galaxies from lower redshift to this epoch. 

Labbe et al. (2009b) recently studied the stellar populations 
of 14 objects at z ~ 7 (all in common with our sample). These 
authors fit models to the stacked WFC3 and IRAC fluxes, de- 
riving a best-fit age of 300 Myr and a best-fit mass of 
x 10 9 M Q , consistent with our results. However, at these 
redshifts nebular [O III] emission could contribute to the flux 
in the IRAC bands, enhancing the amplitude of the perceived 
Balmer break between the observed WFC3 and IRAC fluxes. 

Figure 12 compares our observations to a number of stud- 
ies from the literature, from 2 < z < 7 (Shapley et al. 2005; 
Reddy et al. 2006; Stark et al. 2009, all of which assumed 
a Salpeter IMF), and shows the evolution of the mass of an 
L* galaxy with redshift. For Shapley et al. (2005) and Reddy 
et al. (2006), we computed M1500 for each object, where M1500 
is the absolute magnitude of an object at rest-frame 1500 A, 
and then computed the average mass in bins at M1500 = -19.5, 
-20.5 and -21.5 mag, as well as the 80% range in each bin. 
We then interpolated these values to find the mass at L* using 
the luminosity function of Reddy & Steidel (2009) for these 
redshifts (we split the Reddy et al. sample into two bins at 
z ~ 2 and z ~ 3). Stark et al. (2009) report the mean mass in 
bins of similar absolute magnitude, as well as the 80% range, 
for LBGs at z = 4, 5 and 6. We performed similar calcula- 
tions to derive the mass at L*, using the luminosity functions 
of Bouwens et al. (2007) at these redshifts. We interpolate 
our observations here to the mass at L* using the luminosity 
function of Oesch et al. (2009). 

Comparing our objects to the literature in this consistent 
manner (i.e., at L*, using the values of M* from Oesch et al. 
(2009) and Bouwens et al. (2009c)) allows us to assess how 
the stellar masses of "characteristic" rest-frame UV-selected 
galaxies evolve with redshift. From 2 < z < 6, the mass of a 
L* LBG drops by nearly an order of magnitude, from 10 10 to 
~ 10 9 Mq. Our observations at z ~ 7-8 confirm this trend, 
finding that the masses of L* galaxies at z > 6.5 are < 10 9 

M Q . 

Also shown in Figure 12 are the known stellar masses of 
galaxies selected on the basis of their Lya emission (Lya 
emitting galaxies, or LAEs) taken from the literature, along 
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FIG. 1 1 . — Left: The black curve shows the probability distribution of the stellar masses averaged over all 23 objects with 6.3 < z„hot < 7.5 in our sample. In 
blue is the distribution of best-fit masses (with the numbers labeled on the right vertical axis). All objects in this sample are best-fit with M < 10 10 Mq, with an 
interquartile range of ~ 4 - 9 X 10 s Mq, which is significantly less massive than L* LBGs at z, ~ 3. The average mass distribution confirms this, showing an 
expected mass range from ~ 7 < log(M/MQ) < 10. The red histogram shows the distribution of maximally allowable masses from the two-burst models, which 
peak at the same mass, but extend to log(M/MQ ) ~ 10.0. Thus, even if these galaxies formed 90% of their mass at z ~ 20, they would still be significantly less 
massive than L* galaxies at z ~ 2-3, which have M (at L*) ~ 10 10 Mq (Shapley et al. 2005; Reddy et al. 2006). Right: Same as the left, for galaxies with 7.5 < 
Zphot < 8.6. We reach similar conclusions that these galaxies tend to be lower mass than z ~ 3 LBGs. There is some evidence for evolution from z ~ 7 to z ~ 8, 
as the latter are more likely to have M ~ 10 7 Mq. 

with their combined interquartile range. The relationship be- 
tween LAEs and LBGs is not yet clear, yet it may be likely 
that LAEs form the low-mass and low-luminosity end of the 
LBG population (e.g., Pentericci et al. 2009). However, the 
masses of LAEs do not evolve much with time, evidence that 
they may be progenitors of larger galaxies at subsequent red- 
shifts (e.g., Gawiser et al. 2006; Pirzkal et al. 2007; Finkel- 
stein et al. 2009; Malhotra et al. 2010, in prep). Comparing 
our results to previous studies of LBGs and LAEs, we find 
that the masses of galaxies in our z ~ 7 and 8 samples more 
closely resemble LAEs (or sub-L* LBGs) at any redshift than 
> L* LBGs at z < 6. 

Furthermore, LAEs are also similar to the z ~ 7 and 8 galax- 
ies in that they have blue colors (see § 6). Shimasaku et al. 
(2006) cite evidence that the incidence of strong Lya emis- 
sion in LBGs rises strongly from 3 < z < 6. In § 6 we es- 
timate that the minimum inferred rest-frame Lya EW of the 
z ~ 7-8 galaxies is ~ 70 A, which meets the selection of a 
LAE. At lower redshift, the typically intrinsically more lu- 
minous LBGs tend to be more evolved objects compared to 
LAEs. We do not see many evolved objects at z > 7, thus 
we conclude that the era of more-evolved and massive LBGs 
comes to an end at z > 7. At this early epoch, we are only 
seeing the progenitors later day, more evolved, galaxies. 

6. LYMAN ALPHA EMISSION 

Galaxies selected on the basis of their high Lya emission 
generally correspond to younger, less massive galaxies with 
lower dust attenuation (e.g., Gawiser et al. 2006; Pirzkal et al. 
2007; Finkelstein et al. 2009). Among the z ~ 3 galaxy pop- 
ulation, Shapley et al. (2003) showed roughly 25-30% of the 
predominantly fainter LBG population exhibited Lya in emis- 
sion with sufficient strength to be selected as a LAE via nar- 
rowband surveys (with rest-frame equivalent width > 20 A). 
However, narrowband selected LAEs appear to be different 
than the typically brighter LBGs. LAEs at each redshift have 
on average similar ages and masses, while LBGs are intrin- 
sically more massive, and appear to evolve towards higher 
masses with decreasing redshift. LAEs may become the dom- 
inant population at z > 6, as the incidence of Lya emission in 



LBGs appears to increase with increasing redshift (Shimasaku 
et al. 2006). It is unclear if the increase in Lya emission in- 
cidence among LBGs is due to decreasing dust obscuration 
(e.g., Bouwens et al. 2009a), or if higher-redshift galaxies 
have higher specific star formation rates, and thus would con- 
tain more massive, young stars. In either case, we may expect 
Lya emission to play a prominent role in the SEDs of LBGs 
at high redshifts (e.g., z > 6). 

6.1. Effect on Best-Fits 

Although we have included the effects of Lya emission 
in our model fits, we have no a priori knowledge that these 
galaxies in fact exhibit Lya in emission. However, simply 
allowing a model to have intrinsic Lya emission does not 
ensure that these photons will escape. The amount of Lya 
is strongly model dependent. We have restricted our models 
to constant star formation histories thus there will always be 
some intrinsic amount of Lya. As discussed above, the z ~ 
7-8 galaxies appear to contain little dust, thus any intrinsic 
Lya emission is not strongly internally attenuated. We note 
that the escape fraction of Lya can also depend on the gas as 
well. For example, at such high redshifts, these galaxies may 
be accreting large amounts of gas from the IGM. If the cov- 
ering fraction of the gas around the star-forming regions is 
large, the Lya emission could be scattered into a diffuse halo, 
which would likely not add significant flux to our broadband 
observations. However, recent theoretical work (e.g., Dekel 
et al. 2009; Brooks et al. 2009; Keres et al. 2009) shows that 
gas from the IGM likely accretes along a few discrete fila- 
ments, thus unless we happen to be observing along those fil- 
aments, the star-forming regions would not be obscured, and 
the Lya escape fraction will be high. 

Because our observations do not demand Lya emission, 
we quantify the effect Lya has on our derived model fitting 
parameters by fitting the observed SEDs to models with and 
without emission lines. Figure 13 shows the best-fit age and 
stellar mass results for models with and without Lya emis- 
sion. The models without Lya emission yield younger ages 
and lower stellar masses than when Lya emission is present. 
This is understandable when looking an at object's SED. For 
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FIG. 12. — The stellar masses of L* (characteristic luminosity) LBGs versus redshift taken from the literature are shown as circles. The masses were interpolated 
to the value at L* per the discussion in the text. The remaining symbols are: blue-Reddy et al. (2006); magenta-Shapley et al. (2005); blue, green, yellow-Stark 
et al. (2009) at z = 4, 5 and 6. The estimated mass at L* at z =7—8 from our results are shown by the large orange and red stars. The displayed error bars represent 
the central 80% range of the data, although the small numbers and tight mass ranges of our sample, particularly at z ~ 8, result in small error bars. For that 
reason, we show our individual galaxies as small stars with their uncertainties, which are typically on the order of a factor of ten. At z < 6, there was a slight 
decrease in mass with redshift, which is now extended by our new results. The background gray circles denote stellar masses of Lya emitting galaxies at 3.1 < 
z < 6.5 (Chary et al. (2005), Gawiser et al. (2006), Pirzkal et al. (2007), Nilsson et al. (2007b), Lai et al. (2007), Finkelstein et al. (2007), Lai et al. (2008) and 
Finkelstein et al. (2009)). The gray hatched region denotes the interquartile range of the LAE masses. The masses of the z > 6.3 LBGs studied here are more 
similar to those of LAEs at all redshifts than LBGs at any redshift < 6. 



example, at z ~ 6.5, an object would have Lya in the F105W- 
band. Lya contributes to the amplitude of the color, mak- 
ing it bluer. Removing this strong emission, the model is 
forced to decrease (make bluer) the Y\os - /125 color, produc- 
ing lower best-fit ages, as shown in the left-hand panel of 
Figure 13. Since a younger stellar population is more lumi- 
nous, the model thus requires less mass to fit the observations, 
as seen in the right-hand panel of Figure 13, showing that the 
masses are also less in nearly all cases without Lya emis- 
sion. Including emission lines in the models are thus neces- 
sary when dealing with star-forming galaxies with low dust 
contents (see also Schaerer & de Barros 2009). We conclude 
that our fits including Lya emission are likely to represent 
more closely the physical properties of our sample, and that 
these results are more conservative as they result in higher 
masses and older ages. 

6.2. Implications for Future Surveys 

Many ground-based searches for z > 7 galaxies focus on 
discovering them via the Lya emission line. The image depths 
necessary to detect the observed NIR continuum light are un- 
feasible with most ground-based observatories, but by choos- 
ing to search for Lya emission at redshifts between strong at- 
mospheric emission lines, a number of studies are underway 
which are expected to yield important results (e.g., Nilsson 
et al. 2007a; Hibon et al. 2009; Sobral et al. 2009, Tilvi et al. 
2009 in prep). 

Figure 14 shows the distribution of inferred Lya line fluxes 



from all 35 candidate galaxies in our sample as derived from 
the best-fit stellar population models in §5. As discussed 
above, we expect these objects to have strong Lya emission 
as they appear to be young and unextincted. The expected 
Lya line flux for a typical object in our z ~ 7-8 samples is 
~ 1-2 x 10~ 18 erg s" 1 cm" 2 , although 4 objects have inferred 
Lya line fluxes brighter than 4 x 10~ 18 erg s" 1 cm" 2 . We find 
that 8/31 objects can have fi ya ranging from 4-15 x 10~ 18 erg 
s" 1 cm" 2 when considering our photometric redshift uncer- 
tainties. In contrast, Hibon et al. (2009) reach a limiting depth 
of 8.3 x 10~ 18 erg s" 1 cm" 2 . The ELVIS survey is expected 
to achieve sensitivities of 4.0 x 10~ 18 erg s" 1 cm" 2 (Nilsson 
et al. 2007a). These current and future surveys would not be 
expected to detect many objects similar to those in our sam- 
ple, which would require a ~ four-fold increase in limiting 
sensitivity. Nonetheless, our sample was selected from deep 
WFC3 observations of the HUDF, thus narrowband surveys 
searching a wider area will likely have a much greater chance 
of discovering brighter and rarer objects. 

7. JOINT CONSTRAINTS ON COSMIC REIONIZATION AND 
GALAXY PROPERTIES 

The presence of galaxies at z ~ 7 and z ~ 8 allows us to test 
how these galaxies contribute to cosmic reionization, which is 
believed to be in progress during this time (e.g., Becker et al. 
2001; Fan et al. 2006; Dunkley et al. 2009). To estimate the 
contribution of the galaxies in our z ~ 7 and z ~ 8 samples to 
the budget of hydrogen ionizing photons, we compute the spe- 



22 



Finkelstein et al. 



1000 



a 



o 

X 

<D 
< 



100 




° 100.0 



>- 



o 

X 
<D 

00 
D 



10.0 



1 .0 



0.1 



10 100 
Age incl. Lya (Myr) 



1000 



0.1 



1 .0 



10.0 



100.0 



Mass incl. Lya (10 8 M ) 



FIG. 13. — A comparison of model fitting results with and without Lya emission in the models. The left panel shows the best-fit ages. Without Lya emission 
the models tend to fit the objects with younger ages in most cases. This is understandable, as without a Lya emission line, the model has to fit the entirety of the 
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error bars are excluded for clarity. 



cific luminosity density of each object at rest-frame 1500 A. 
For each object, we interpolate the WFC3 T105, /125, and H\60 
photometry and apply the luminosity distance to derive the 
specific luminosity at 1500 A. We then sum the individual lu- 
minosities of each galaxy and divide by the comoving volume 
element over 6.3 < z < 7.5 and 7.5 < z < 8.6 for the z ~ 7 
and z ~ 8 samples, respectively. We make no corrections for 
survey incompleteness or for galaxies fainter than the magni- 
tude limit of our catalog, although we note that either of these 
corrections would increase the luminosity densities we cal- 
culate. Therefore, our luminosity densities are conservative 
lower limits. 

Figure 15 shows the rest-frame 1500 A specific luminos- 
ity density we derive for our samples at z ~ 7 and z ~ 8. 
We estimate errors on the luminosity densities by performing 
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higher, with 8 objects having i Lya > 4 X 10 18 • 



a bootstrap simulation whereby we generate galaxy samples 
by randomly selecting galaxies from the true sample with re- 
placements including the effects of Poisson uncertainties and 
the photometric uncertainties. 

We compare the luminosity densities we derive to models 
for the required number of ionizing photons to keep the Uni- 
verse ionized at a given redshift (Madau et al. 1999). Follow- 
ing Madau et al. (1999) and Pawlik et al. (2009), we find the 
required luminosity density at 1500 A needed to maintain an 
ionized universe is 



puv = 1.25 x 10 25 e 53 ( ^P- 



0.0463 



C 

Js 



-z-, (5) 



in units of erg s 1 Hz 1 Mpc 3 , where f2& is the cosmic baryon 
density and h-jo is the Hubble parameter in units of h = 0.7, 
and the constant 0.0463 is our assumed value of fib (x /i 2 ) 
from Wilkinson Microwave Anisotropy Probe Year 5 data (Ko- 
matsu et al. 2009). The variable £53 is the number of Lyman 
continuum photons per unit of forming stellar mass is units 
of 10 53 photons s" 1 (M yr" 1 )" 1 . To determine this value, we 
used CB07 models with a Salpeter IMF, a constant SFH and 
Z = 0.2 Zq, which is within the 95% confidence range on the 
metallicities of our objects from §5, and gives £53 = 1.2. Vary- 
ing the metallicity fromZ Q - 0.02 Z changes £53 to 0.9 - 1.4, 
respectively. We note that if the high-mass end of the IMF 
is biased towards more massive stars, more photons would be 
emitted at 900 A relative to this equation, thus our assump- 
tions are conservative. 

The equation above depends on two variables, the first 
being the "clumping" factor of neutral hydrogen, C = 
(Phi) (Phi) -2 - The second variable / esc is the average escape 
fraction of hydrogen-ionizing photons. While Eqn. 5 uses 
a single value for f esc this is the average over all galaxies, 
and certainly there is a distribution, possibly with a luminos- 
ity dependence. Neither C nor f esc are well understood theo- 
retically or empirically. While Madau et al. (1999) argued for 
a high H I clumping factor C w 30, recent simulations and ra- 
diative transfer calculations imply a lower average clumping 
of C « 5-6 (Pawlik et al. 2009; Finlator et al. 2009), due to 
photo-heating of the IGM by star formation. The escape frac- 
tion of ionizing photons has been measured at lower redshifts, 
with / esc < 0.1 at z = 0, /esc = 0.08 at z ~ 1 (Siana et al. 2007), 
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FIG . 15. — The circles represent the specific luminosity density of our sam- 
ple, uncorrected for dust or incompleteness, split into two bins in redshift. 
We denote the regions where we believe the Universe is mostly ionized and 
mostly neutral. We are now beginning to probe into the region where reion- 
ization is likely to be beginning. The curves show the critical specific lu- 
minosity density needed to reionize the Universe for ratios of the clumping 
factor C over the escape fraction f esc . If our points lie above a given curve, 
that indicates that the Universe would be fully reionized by that redshift. The 
downward triangles denote the total specific luminosity density for our sam- 
ples if we account for fainter galaxies down to M1500 = -15 with the luminos- 
ity function of Oesch et al. (2009). For large values of / E jc(~ 0.5) galaxies 
at these redshifts are able to reioinize the universe. Given the evidence for 
low dust extinction, especially from the fainter galaxies in our sample, high 
values of f esc may be likely. 

and measurements of / esc = 0.15-0.6 at z ~ 3 (Steidel et al. 
2001; Shapley et al. 2006). The escape fraction at z > 6 is 
unknown, but current models require f esc > 0.2 (Pawlik et al. 
2009; Finlator et al. 2009). 

Figure 15 shows several curves from equation 5 for various 
combinations of C/f esc . Based on this figure we conclude that 
galaxies at z ~ 7 and z ~ 8 have sufficient hydrogen-ionizing 
photons to maintain cosmic reionization under certain condi- 
tions. Even with no corrections to our measured luminosity 
densities, the galaxies in our sample would reionize the Uni- 
verse for C//e. sc = 1 to ps 10. This ratio of C to / esc already 
encompasses theoretically motivated values. For example, 
taking C = 5 from the recent radiative-transfer calculations 
of Finlator et al. (2009) implies that for / esc = 0.5, the galax- 
ies we observe would nearly reionize their local volume by 
z~6.5-7. 

Although we have considered only the UV luminosities of 
the galaxies in our sample, we can estimate the contribution 
from fainter galaxies using the z ~ 7 UV luminosity function 
of Oesch et al. (2009) with our z ~ 7 sample, and shifting it 
to match our z ~ 8 sample. We calculate a correction factor 
as the ratio between the integral of the Oesch et al. luminosity 
function down to M1500 = —15 by the integral of our observed 
luminosity function at z ~ 7 and 8. The implied corrections 
are a factor of 3 at z ~ 7 and 4 at z ~ 8. These are indicated by 
the downward triangles in Figure 15. These points represent 
plausible upper limits on the luminosity density (barring addi- 
tional corrections from incompleteness and effective volume). 
Achieving higher luminosity densities from these galaxies 
would require substantial corrections from galaxies fainter 
than our magnitude limits (c.f. Yan et al. 2009). Allowing 



for the maximal hydrogen-ionizing photon escape fraction, 
fesc = 1 1 the dumpiness must be C < 20 to match the measured 
UV luminosity density and still maintain a partially ionized 
IGM. For more plausible escape fractions, /esc ~ 0.1-0.5, 
the clumping factor must be closer to C =3-5. 

The low levels of dust extinction inferred from the very blue 
rest-frame UV colors of our sample implies that the escape 
fraction for hydrogen-ionizing photons may be quite large 
(though see §6.1). Indeed, although our above calculations 
have assumed no dependence of f esc on the UV luminosity 
(and thus the stellar mass as the stellar mass correlates with 
UV luminosity at these redshifts - see Tables 1 and 3), this 
may be the case, as we showed above that the fainter galaxies 
may be slightly bluer. Therefore, it may be that the fainter 
galaxies have higher ionizing escape fractions, / e sc > 0.5, re- 
sulting from combination of low metallicity and minimal ex- 
tinction in these galaxies. As galaxies develop higher metal- 
licity and dust content, it may be they achieve lower escape 
fractions, consistent with measurements at z = 3 (Shapley et al. 
2006), but by then the universe is well ionized. 

7.1. Comparison to Recent Studies 

Several recent studies have also made estimates of the con- 
tribution of galaxies at z ~ 7- 10 to the UV luminosity density 
required to reionize the Universe. Labbe et al. (2009b) used 
average WFC3 and IRAC colors from stacked data for their 
sample to estimate the average SFR and stellar mass. Assum- 
ing the stellar mass in their galaxies formed at 7 < z < 11, 
they arrive at a result consistent with the one we report here, 
assuming similar values for the H I dumpiness and UV escape 
fraction. Similar results are achieved by Bunker et al. (2009) 
and Yan et al. (2009), if we adopt our values of C and / e sc (and 
we note that the extremely steep slope of the UV luminosity 
function derived by Yan et al. is able to reionize completely 
the Universe at z ^ 7). Therefore, the following (general) 
conclusion is fairly robust, as it has been derived from several 
studies using different techniques. The galaxies we identify 
at z ~ 7 and z ~ 8 produce a sufficient UV luminosity den- 
sity for reionization. The specific details are still unclear (in 
particular the unknown values for the UV escape fractions), 
and will require more observational studies of high-redshift 
galaxies. Based on our study, and those in the literature, we 
predict that the galaxies at z ~ 7 - 8 have / esc ps 0.1 — 0.5 in 
order for reionization to proceed. 

8. CONCLUSIONS 

We have analyzed new, very deep data from WFC3 on 
board HST of the Hubble Ultra Deep Field, obtaining a sam- 
ple of 31 galaxies at 6.3 < z p hot < 8.6. We have examined the 
colors, physical properties, and Lya-emission characteristics 
of these galaxies using a slate of empirical galaxy templates. 
We find that the rest-frame UV colors of these galaxies are 
very blue, bluer than the expected colors of local starburst 
galaxies (> 4a significance). This represents strong evolution 
from z ~ 3, where typical star-forming galaxies have colors 
consistent with these local starbursts. Taking into account a 
full simulation of the galaxy photometric uncertainties, we 
find little evidence that the galaxies in our z ~ 7 and 8 sam- 
ples host exotic stellar populations with primitive metallici- 
ties. Rather these objects appear similar to very blue (nearly 
unextincted) local starburst galaxies such as NGC 1705. 

We fit the measured colors of the galaxies in the z ~ 7 and 8 
samples to stellar population models, using the updated code 
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of Chariot & Bruzual (2007). We find that although the con- 
straints on the ages and dust extinction are marginal, the data 
are consistent with the stellar populations in these galaxies 
having low ages (^100 Myr), little-to-no dust extinction, and 
low metallicity (< 0.05 Z Q at ~ 1 a confidence). 

We derive tighter constraints on the objects' stellar masses, 
ranging from 10 s — 10 9 Mq for z ~ 7 and 8 galaxies. These 
results extend the evolutionary trend that the stellar mass of L* 
galaxies decreases with increasing redshift. We examine the 
maximal allowable stellar masses in these objects by fitting 
them to a two-burst model, which includes an old "burst" of 
star formation with 90% of the total stellar mass, and with an 
age equal to the time from z = 20 to the photometric redshift 
of each galaxy. The young age of the universe at these high 
redshifts (look back time of t = 750 Myr at z = 7) places strong 
constraints on the amount of old stars in these objects, much 
more so than the upper limits derived from the deep IRAC 
data. With this age constraint, even the maximal masses for 
the z ~ 7 and 8 galaxies are typically only a few times 10 9 M Q . 
These are larger than their best-fit masses using single stellar 
populations, but they are significantly less than those derived 
for "L*" star-forming galaxies at z < 6. 

We investigate the effect of Lya emission on the interpre- 
tation of the colors of the galaxies in the z ~ 7 and 8 sam- 
ples. As may be expected for these relatively unextincted, 
star-forming objects, Lya can be strong. If this emission is 
not included in the models, the best-fit results systematically 
underestimate both the age and stellar mass, as the models 
become bluer to account for the observed rest-frame UV col- 
ors. We estimate the Lya emission in these objects from the 
best-fit stellar population models, finding that while all ob- 
jects have rest-frame EW >70 A, the implied line fluxes for 
typical z ~ 7 and 8 galaxies are below the sensitivity of cur- 
rent and planned Lya searches by factors of a few at these 
redshifts. 

Additionally, the derived stellar masses for the z ~ 7 and 
8 galaxies are comparable to those of galaxies at 3 < z < 7 
selected on the basis of their Lya emission, which tend to be 



less evolved than typical lower-redshift LBGs. We conclude 
that evolved galaxies at z > 7 are very rare, and that the galax- 
ies in our samples represent the progenitors of lower redshift 
galaxies in their first throes of star formation. 

We calculate the UV luminosity density of the galaxies in 
ourz ~ 7 and 8 galaxies, and compare it to the necessary value 
for cosmic reionization. We find that the galaxies in our sam- 
ple provide sufficient UV luminosity density (within factors of 
a few) of that required for reionization, even without making 
corrections for incompleteness. Estimating the contribution 
from galaxies below our detection limit, star-forming galax- 
ies at z ~ 7 and 8 can ionize the the IGM if they have UV 
escape fractions / esc w 0.3-0.5. These high escape fractions 
may be likely given the the low dust extinction and metal con- 
tent in these galaxies. 

The recent installment of WFC3 on board HST now allows 
for the study of faint objects (~ 29 mag) in deep near-IR 
imaging. Even with the small amount of data available for 
these galaxies we are able to derive useful constraints on their 
properties, and their effect on the IGM. Future infrared data 
will allow better constraints on the physical properties of these 
objects to link them both to lower redshift galaxies, and to 
galaxies yet to be discovered at higher redshifts. 
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